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The Infrared Absorption Spectra of Nitrogen Dioxide and Tetroxide* 
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The absorption spectra from 1 to 4u of pure NO.-N2O, 
mixtures in a long tube were examined with a rocksalt anda 
quartz spectrometer. Several new bands were discovered 
from 1.9 to 4y in the spectra of both NOz and N2O,. The 
NO: bands were analyzed as combinations of two funda- 
mentals. The N2O, bands were interpreted as combinations 
and harmonics of reported fundamentals. Three bands of 


N20, lie at shorter wave-lengths than correspond to the 
work of dissociation into 2 NO2. A study of the symmetry 
properties shows that only those models where the NO2 
groups are joined by the nitrogen atoms would give rise to 
the observed infrared spectrum. Some features of the 
spectrum seem to be in favor of the existence of molecules 
where the NO» groups do not lie in the same plane. 





INTRODUCTION 


HE infrared absorption spectrum of the 
NO.-N2O, system has recently been in- 
vestigated by several authors.!: *: 3» 4 We present 
in this paper the absorption spectrum, in the 
nearer infrared, of deep layers of some very 
pure NO», and N2O,,° as resolved by a rocksalt 
and by a quartz prism spectrometer. The spec- 
trum of NO, in this region is composed of distinct 
bands which we consider to be combinations of 
the fundamentals reported by Sutherland.!: 6 
Two bands of NO: are discovered which we have 
interpreted as combinations of some undiscov- 
ered fundamental with one already known. 


* Contribution No. 319. 

1G. B. B. M. Sutherland, Proc. Roy. Soc. A141, 342 
(1933). 

* R, Schaffert, J. Chem. Phys. 1, 507 (1933). 

7A. B. D. Cassie and R. L. Bailey, Nature 131, 239 
(1933). 

* J. Strong and T. Woo, Phys. Rev. 42, 267 (1932). 

*The preparation will be described elsewhere by L. 
Harris and W. S. Benedict. 

*°G. B. B. M. Sutherland, Proc. Roy. Soc. A141, 535 
(1933). 


APPARATUS 


Since the gases react with cements, metals and 
greases, they were confined in an all quartz ab- 
sorption tube 100 cm long and 2.54 cm diameter, 
surrounded by an electric furnace. The windows 
were of carefully selected quartz and, although 
only 2.0 mm thick, were opaque beyond 4u. 

In order to reach this wave-length no other 
quartz was permitted in the optical path, in one 
series of experiments. The light of a Nernst glower 
was reflected by gold mirrors through the tube 
and then through a minimum deviation 60° rock- 
salt prism spectrometer. A Coblentz twenty-two 
junction thermopile traversed the spectrum with 
twice the angular velocity of the prism and the 
angle of (minimum) deviation of the incident 
light was measured by a graduated table reading 
to five seconds. This spectrometer was calibrated 
by the simple minimum deviation formula. 
Paschen’s’ values for the refractive index of rock- 


_ salt corrected to 21°C were used but so few points 


51 


are given in this region of inflection that a con- 


siderable uncertainty in the wave-length arises 


7 Paschen, Ann. d. Physik [4] 26, 120, 1029 (1908). 
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(+0.003u). The prism was kept at constant 
temperature within 0.1°C since dd/dt=0.015y/ 
degree. 

In the region of transparency of quartz 
(<2.8u) much higher dispersion was obtained by 
means of a quartz half-prism constant deviation 
spectrometer of large aperture. The light of a 
1000 watt tungsten lamp was passed through the 
meter tube by quartz lenses. Monochromatic 
light from the exit slit was focussed on a Kipp 
thermopile. A graduated drum measured the 
angle between the half prisms in radians but the 
constants of the calibration had to be evaluated 
by locating mercury emission lines (to 1.71,). 
Calibration from the positions of the bands of 
common substances did not seem advisable as the 
values in the literature do not agree. 

With each instrument the positions of the ab- 
sorption bands of water-vapor® in the air and of 
three NO: and N.O, bands? agreed well within the 
experimental error of 0.0124 with the published 
values. 

The intensity of the transmitted light was 
measured by the thermopiles and a Leeds and 
Northrup HS galvanometer, giving reproducible 
transmission curves within one percent. Two to 
ten readings were taken per slit width (see 
figures) and the curves determined several times 
at each absorption band. The tube was emptied 
by freezing out the contents into a side arm with 
“dry ice’’ (commercial solid CO2). The transmis- 
sion of the full tube was found at 21°C and 130°C 
with the rocksalt instrument and at 24°C and 
160°C with the quartz, where the calculated com- 
positions were as follows: 

Temperature of ab- 
sorption tube 

Temperature of res- 
ervoir 

Partial pressure re- 


duced to 0°C 
of pNO2(Pxoy) mm 


g 
of Ns0«(Px0,) 
mm Hg 


120 130 160 
—12 50 75 


(°C) 21 24 
(°C) 21 24 

188 217 188 1034 
511 495 2 86 


TREATMENT OF DATA 


Although we are not primarily interested in the 
magnitudes of the absorption coefficients, it is 


8G. Hettner, Ann. d. Physik [4] 55, 476 (1918). 
9 See tables below. 
10 M. Bodenstein, Zeits. f. physik. Chemie 100, 68 (1922). 


necessary to resolve the absorption curve for each 
equilibrium mixture into its two components in 
order to locate the bands of each constituent. Us- 
ing logarithms of transmissions, a correction was 
applied throughout the spectrum to the absorp- 
tion curves (rocksalt data) of the mixtures by 
subtracting the absorption due to the minor 
constituent, calculated by Beer’s law from the 
absorption curve of the other mixture where this 
species was now in excess. It was assumed that the 
absorption coefficients do not depend on the 
temperature or pressure under this resolution; 
later experiments with NO, showed Beer’s law 
was valid within our experimental error, as 
might be expected from the fact! that the funda- 
mentals even under high resolution show no de- 
tailed structure. 

With the quartz spectrometer, the absorption 
tube was heated to 160°C, where the gas was 
practically pure NO», so that no correction was 
applied for the absorption of the small amount of 
N20, present. To find the absorption of pure N.O, 
we had first to estimate the absorption of the NO, 
in the gas mixture at 24°C by another series of 
measurements made with the tube at 120°C and 
the reservoir at — 12° where the gas was all NOz, 
and at the same pressure as in the mixture 
at 24°C. 

Fig. 1 shows the curves for NO, at 130°C and 
N2O, at 21°C as obtained with the rocksalt 
spectrometer. Fig. 2 gives the absorption of NO, 
at 160°C and NO, at 24°C under the higher dis- 
persion, all curves having been corrected as out- 
lined above. A comparison shows how the ap- 
parent absorption increases with the resolution 
and how overlapping bands shift the maxima. 


RESULTS 

The absorption coefficients of NOz from 1 to 
2.1 and of N2O, from 0.75 to 1.9 are less than 
0.0003 atm.—'cm™ (three percent being the ac- 
curacy of measurements of weak absorption), 
unless there are bands narrower than the slit 
width (shown in figures). The spectra of both 
gases between 1.9 and 4y consist of more or less 
isolated bands whose maxima are given in Tables 
I and II. 

We use Mecke’s notation, which is identified 
with the notation Sutherland used for the funda- 
mentals in a later table. If the z-axis be taken 
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TABLE I. NO» at 130° and 160°C. 








Ob- 
served 
by others 
cm! 


648!: 2.3 mo ox Strong 
(1321)5 1321 ve — 
1621!-2-3 1621 wy, Very strong 
2220! 2236 =vg+dr Medium 


26305 2630 2» Medium 

2910! 2893 vette Very strong 

3242! 3206 92, No maximum 
but strong 
absorption 


Vo +r +6, Weak 
3vx Hidden weak 


maximum 
4140 Veo t2vez 


Strong 
4560 Qvetve Hidden weak 
4753 3u¢ 


Calcu- 
lated 
cm7! 


Analysis Intensity 





{ 2601 
2667 

2905 
ca. 3200 


3597 


{ 3454 
ca. 3930 


maximum 


Medium 








TABLE II. N20, at 21° and 24°C. 








Ob- 
served 
by 
others 
cm7! 


3804 ca 
682! bo ‘a 
752' ona 
8804 bea +5es 
ol Voxs +6, ‘a 
1 Vra 
1744} res 
1750 Vra +5es 


Vea tors 
2060! Vrs(2) +bra 


VastOra appa 
Voat brs Medium 


Strong 
Vara tvrs(2) Strong 
Vos +Vra 


perp. ve tVra 
3077 


3091 Vea +yzs(2) 
3478 Voa +e 
perp. 2v, 
See 

2ves +Vra 

YVoa +2vz0 

Vos +vratvrs 
Voa +2ves 


2vee +¥x0 

2¥¢a +7r0 
perp. 2vg+vra 
VoatVestVrs 
perp. 2vg+vrs 
2v¢s +¥ea 

3¥gq OF 3r¢ 


Analysis Intensity 





2501 
{ 2608 
2649 
2967 
Very strong 


Strong 


Weak 
Strong 
Medium 
Weak 
Weak 
Weak 


Medium 


3861 
4160 
4230 
4340 
4478 
4620 


4706 
4778 
5197 


Medium 
Medium 








* Anharmonicity correction (not attempted) would 


justify this assignment. 


KING 


through the N—N bond, then z= signifies || to 2; 
o, | toz;o’, | tozand L tothe plane of the NO, 
groups; the vibrations of the NO: groups will be 
symmetric (s) or antisymmetric (a) with respect 
to the center of symmetry. The assignment is for 
the plane molecule; for the perpendicular model 
the assignment is the same, except that the (os) 
and (ca) vibrations become identical, as will be 
discussed below. 


THE SPECTRUM OF NO, 


Only two fundamentals have been reported for 
NOs, at 1621 and 648 cm!.!: 2: 3 If the molecule 
is linear the third frequency v, can be calculated 
to be about 1000 cm and although the funda- 
mental would be absent in the infrared, combina- 
tions should be expected. The absence of a third 
frequency however is not a definite indication of a 
linear molecule; there are several triangular 
molecules known where the v, frequency appears 
very weakly. Harris and Benedict’ have found 
separations in the ultraviolet spectrum of NO; 
of 1321 and 2630 cm-, which are very reasonable 
values for v,, 2v, of a triangular molecule. These 
values fit the infrared data very much better than 
a calculated 1000 cm~. If », is weak, the other 
overtones, U.V, Ve¥etve, etc., (ve even) will also 
be weak (cf. H,O,"). Table I shows that indeed 
the only strong bands are the o-bands of odd 1,, 
their sharpness warranting the evaluation of 
three anharmonicity terms. The frequencies were 
calculated from 


E= 1681(v, +3) +1357(v,+3) _ 18(v, +3)? 
— 6(v¢-+})?—48(0e-+4) (Ce +4). 


The intensity of »,+v, and v,+2», is reflected in 
the magnitude of the interaction constant, 
—48 cm". 

Bands reported! * at 1350-75 cm are prob- 
ably the frequency 1365 cm! of nitrite ions’ 
formed by reaction of the gases with the rocksalt 
windows, avoided in the present work. These 
bands may have masked the weak NO, funda- 
mental at 1321 cm“. 


THE SPECTRUM OF N,O, 


The new bands between 1 and 4y can be satis- 
factorily interpreted as combinations of the re- 


4 Mecke et al., Zeits. f. Physik 81, 313 et seq. (1933). 
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TABLE III. 


TABLE IV, 








Operators 


E C2(z) C2(y) Co(x) Ci Cs(z) Ce(y) Ca(x)| Configurations 


Operators 


E C2(z) 2C;C4, —4(2) 2C2(x, y) 2Cs’(x’, ’) Configurations 





Zero state 
Yrs=M= eed 
6x3 = vs =813 
veges = vs= 283 
Free rotation v7 


1 |Z, vea=v1=1265 
ra— ¥3> 752 
6q’s=vi2 
Y, vea=vs= 1744 
bea = v9= 380 
Yog=ve=1724 
b¢s = 1c = 500 
X, 66’ = 111 = 682 














ported infrared and Raman fundamental fre- 
quencies. As was the case with the data treated 
by Sutherland, the comparative scarcity of bands 
in this region (if combinations and overtones of 
all twelve fundamentals appeared the absorption 
at the present dispersion would appear nearly 
continuous) is positive evidence for a symmetrical 
N.O, molecule, undoubtedly one where the NO, 
groups are joined by the nitrogen atoms. How- 
ever, we must examine the data for indications of 
two well-populated configurations of N2O, of 
nearly the same energy, in which the NO2 groups 
lie in the same plane or are mutually perpendicu- 
lar, and for indications of free rotation through 
both these positions upon surmounting a com- 
paratively low hill of potential energy. 

Neglecting for the present the interaction with 
free rotation, the other eleven fundamental fre- 
quencies may be arranged among the various 
representations of the symmetry groups corre- 
sponding to the equilibrium positions. The selec- 
tion rules of combination, derived from sym- 
metry considerations,’ '* are equally rigorous 
for the two configurations and are equally com- 
patible with the observed spectra. 

Plane N2O, has the symmetry properties of 
group V; which has the irreducible operators and 
representations given by Table III; the z-axis is 
chosen as the N—N bond; the y-axis lies per- 
pendicular to this in the plane of the molecule. 

The selection rule for the infrared, namely 


f (am /aq.) Winie. (£)W rinar(€)4 


must be positive in all symmetry operations. 


* L. Tisza, Zeits. f. Physik 82, 48 (1933). 
*R.S. Mulliken, Phys. Rev. 43, 279 (1933). 


1 1 Zero state 
Vrs = 1360 
vigg= 283 
bxs= 813 


Free rotation 
¥ra= 1260 
bxra= 752 


x, ¥- 
vg = 1744 
6¢ =500382 
bg’ = 682 














Applied to this group, it shows that for transitions 
from the zero state (1) no fundamentals of synt- 
metrical (s) vibrations, or (2) even overtones of 
any fundamental are active and (3) there are no 
combinations of s with s or of a with a or (4) of 
vibrations lying in representations with the first 
subscript the same. The most intense of the 
permitted combinations have been given in 
Table II. 

The perpendicular molecule is of the Vg group, 
described by Table IV. 

The outstanding difference between the groups 
is seen to be the occurrence in Vz of the degen- 
erate representation Ey, wherein the 6 o-vibra- 
tions coalesce to form three (not two as stated by 
Sutherland) degenerate pairs, which are all ac- 
tive. However, as only the three (ca) are active 
in V,, there will be the same number of observ- 
able infrared fundamentals as in Vg. That the 
values of the Raman frequencies® do not coincide 
with those of the infrared for these pairs cannot 
be taken as evidence against the existence of the 
perpendicular form in the gas, as the Raman 
spectrum was obtained only in the solid, which 
x-ray studies'* have shown to be composed of 
predominantly plane molecules. In the overtones 
and combination bands also, Va gives the same 
number of active bands as Vy. In Va, even- 
numbered harmonics of the o-vibrations will 
occur (with pronounced Q branches) but with 
less than the full intensity expected of states of 
increased (degenerate) weight, since the charac- 
ters of the product, £,*, split into the sum 
A,+B,+B2, of which only the last is active, the 
representation of the Z-moment. All combina- 
tions of the o-vibrations with each other and with 


other types of vibrations are also permitted. 


4S. B. Hendricks, Zeits. f. Physik 70, 699 (1931). 
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Throughout, a correlation of the permitted com- 
binations of Va with V;, may be made. 

In Table II we have shown how the permitted 
combinations, on either model, account for the 
main features of the spectrum. The following 
considerations, however, lead us to suspect that 
both models participate in giving the observed 
bands. 

The point at issue is whether the o-vibrations 
have one frequency or two. Only », is observed 
with certainty in overtones and combinations. In 
the region of 1750 cm~ the envelope of absorp- 
tion as observed by Sutherland is wide and asym- 
metrical with a deep minimum which Sutherland 
interpreted as the central gap to be expected in 
the fundamental »,, where the moment vibrated 
along the middle axis of inertia of an asym- 
metrical top and hence as evidence for the plane 
molecule. A rough calculation with the aid of 
Dennison’s diagrams!: ° reveals, however, that 
for a molecule as heavy as N2O, the central gap 
would not be so wide as observed and that in all 
probability the envelope, under the attainable 
resolution, would appear with a wide and weak 
central maximum (ratio of moments 1 : 1.8 : 3.1). 
In the perpendicular model, N2O; is a sym- 
metrical top (ratio 1 : 2.6 : 2.6) and would havea 
somewhat stronger central maximum. It is there- 
fore probable that the observed absorption is due 
to the fundamental v.,4)=1744 cm, on which 
is superposed the permitted combination y,, 
+6,.s)=1765 cm. The Raman frequency of the 
solid (plane) occurs at 1724 cm~; hence we must 
look for combinations of 1744 and 1724 from the 
plane and for combinations of 1744 only from the 
perpendicular molecules; if rotational oscillation 
is present the bands will lie between these limits. 
In the region of the fundamental the presence of 
weak asymmetric absorption lying to shorter 
frequencies of 1744 may be due to the 1724 (v,.) 
coupling with v;o¢ to become weakly active. 

The first ‘‘overtone”’ would lie at 3468 cm= 
(less some anharmonicity term which must be 
present if the band is to appear with appreciable 
intensity) if plane, and at 3488 cm if perpen- 
dicular. The observed band at 3478 cm is pre- 
sumably the center of the overlapping of a plane 


% Dennison, Rev. Mod. Phys. 3, 280 (1931). 
16 Gerhard and Dennison, Phys. Rev. 43, 197 (1933). 
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model combination, vea+ vs. = 3468 cm with the 
perpendicular frequency 2v, = 3488 cm—. 

The calculated separation between 2y¢(a)+vr0 
and e(ay+¥o(s)+¥ss is 75 cm for the plane, 
95 cm~ for the perpendicular molecule. The ob- 
served is 70+5 cm but a correction for slit- 
width would increase this value. Moreover the 
envelope at 4778 cm™ seems to be composed of 
the bands of both types of molecules. 

All these facts are evidence for the existence of 
an appreciable fraction of molecules in the gas in 
which the NOs» groups do not lie in the same plane. 
It would be of interest to examine the envelopes 
of these overtone bands under high resolution in 
the attempt to decide whether the plane or the 
perpendicular form is the equilibrium position 
of the gaseous N.O, molecule and to determine 
the height of the potential hill hindering free ro- 
tation. However, the discussion we have given 
should make it clear that an unambiguous deci- 
sion will be attended with difficulties. 

As the single N—N bond should have no 
directional properties and the forces between the 
rather distant O atoms on opposite NOz groups 
are doubtless small, a considerable fraction of the 
molecules, even at low temperatures, will prob- 
ably be oscillating widely about the equilibrium 
position, if not rotating freely. The oscillation 
motion will be strongly coupled to the other 
modes of vibration in such a way that these 
deviate from the frequencies and selection rules 
typical of the equilibrium position towards the 
frequencies and selection rules of the less sym- 
metrical configurations taken up during the 
oscillation. For completely free rotation, the 
coupling will be so strong that the other vibra- 
tions will resemble exactly neither those of the 
plane nor the perpendicular models but will be 
intermediate between them. The long stretches of 
continuous absorption observed in certain regions 
may indicate some degree of free rotation. 


DISSOCIATION OF N2O, 

The energy of dissociation of NsO, into 2NO: 
as given by equilibrium constants” is 12.6 kcal.. 
equivalent to 4440 cm-. N,Q, absorbs light of 
shorter wave-length than corresponds to its en- 
ergy of dissociation in three strong, well-defined 
isolated bands, centers at 5197, 4778 and 4706 
on™. 
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ABSORPTION SPECTRA OF NO, 


We have analyzed these bands as combinations 
of vibrations within the NO, groups, correspond- 
ing to the bands 3», and 2»,+ », of single NO. and 
not of the vibration of the N—N bond (v,,). Asa 
result of this absorption the molecule is elevated 
to a quantized level of a potential surface above 
the convergence limit of the decomposition fre- 
quency (v2x,;). The molecule may now either re- 
distribute its energy into the N—N bond and de- 
compose by a radiationless transfer through the 
col above the convergence limit or it may, more 
probably, lose this energy by collisions. The ac- 
tual occurrence of discontinuous absorption is of 
interest from the standpoint of the old ‘“‘radiation 
hypothesis” and as a verification of a phenom- 
enon that has long been suspected to be possible. 

Even if all the energy absorbed in this region 
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resulted in the decomposition of the molecule, the 
amount of infrared light needed to oppose meas- 
urably the thermal rate of formation is exceed- 
ingly high at ordinary temperatures (10® mols/sec. 
requiring 10° kcal./sec.). It may however be pos- 
sible to decompose, at low temperatures, the 
white vitreous solid N,O, into red NO, by infra- 
red radiation between 1.92 and 2.25y. 
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Crystal Structures and Expansion Anomalies of MnO, MnS, FeO, Fe;0, Between 
100°K and 200°K* 


BENNETT S. ELLEFSON{ AND NELSON W. TAyLor,{ School of Chemistry, University of Minnesota, Minneapolis, Minnesota 
(Received November 18, 1933) 


(1) Determinations have been made of structure type 
and lattice constants of MnO, MnS, FeO and Fe;Q, at 
temperatures near those at which these substances show 
specific heat anomalies. (2) No major structure changes 
were apparent in the x-ray patterns. All structures are 
cubic at all investigated temperatures. Lattice constant 


data indicate similarity in the abnormal expansion behavior 
of MnS and FeO and of MnO and Fe;QO.. (3) The specific 
heat and expansion anomalies are probably due to electron 
transitions within certain atoms in the crystal lattice and 
may be analogous to the behavior of metallic iron at 
765°C, or of quartz at 575°C. 





INTRODUCTION 


CCORDING to results of recent experi- 

ments carried out by Millar! and Anderson? 
at the Pacific Experiment Station of the Bureau 
of Mines, MnO, FeO, Fe;0; and MnS show 
anomalies in their specific heat curves. Each 
curve shows a cusp or hump superimposed upon 
the normal ascending trend. The plots of heat 
capacity versus temperature for these four com- 
pounds are reproduced in Figs. 2, 3, 4 and 5. 
For exact numerical data reference should be 
made to the original publication. 

Whenever a transition in state occurs, there 
results an abnormal absorption or evolution of 
energy. Consequently, the relation between heat 
capacities and temperature is quite helpful in 
locating such transitions. The change from one 
structure to another is usually accompanied by a 
heat absorption or evolution over a relatively 
short temperature range. 

Experiments have been reported which show 
that the abnormal absorption of energy is not 
always accompanied by change in crystal struc- 


* Presented by Nelson W. Taylor before the Division of 
Physical and Inorganic Chemistry at the Washington 
meeting of the American Chemical Society, 1933. 

t Submitted in partial fulfillment of the requirements 
for the degree of Master of Science, June 1933. 

t Present address: Head, Department of Ceramics, The 
Pennsylvania State College, State College, Pa. 

1R. Millar, J. Am. Chem. Soc. 50, 1875-83 (1928); 51, 
215-22 (1929). 

2C. T. Anderson, J. Am. Chem. Soc. 53, 476-83 (1931). 


ture. Simon* discovered anomalous specific heat 
behavior in ammonium chloride at some point 
between 223° and 243°K. Investigation of the 
crystal structure at these low temperatures by 
Simon and v. Simson‘: > shows anomalous be- 
havior at 242.3°K for ammonium fluoride, at 
242.6°K in the case of ammonium chloride, at 
235°K in the case of ammonium bromide and at 
230.5°K in the case of ammonium iodide. 
Klinkhardt® in his work on specific heat of iron 
found two outstanding abnormalities in the 
specific heat curve. One of these abnormalities, 
occurring in the neighborhood of 906°C, repre- 
sents a sharp drop in the specific heat curve and 
is explained by Westgren and Lind’s’ work 
showing a change from the body-centered cubic 
lattice of the alpha iron to the face-centered 
cubic lattice of the gamma iron. The other 
specific heat anomaly which takes place over a 
considerable temperature range and_ which 
reaches a maximum at 765°C, cannot be ac- 
counted for on the basis of crystal structure 
change. However, it is interesting to note that 
the peak of this curve occurs at the same point at 
which the iron loses its ferromagnetic properties.” 


’F,. Simon, Ann. d. Physik 68, 241 (1922). 

4F, Simon and v. Simson, Naturwiss. 38, 880 (1926). 

5 F. Simon and v. Simson, Zeits. f. physik. Chemie 129, 
339 (1927). 

6H. Klinkhardt, Ann. d. Physik 84, 167-200 (1927). 

7 A. Westgren and A. Lind, Zeits. f. physik. Chemie %. 
181 (1921). 

8 Weiss, J. de Physique 7, 249 (1908). 
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CRYSTAL 


APPARATUS 


The source of x-rays for all pictures taken in 
this experiment is a Siegbahn metal x-ray tube 
with an iron target. The camera is one designed 
by Taylor.’® Slight modifications were made to 
reduce thermal capacity. The specimen is 
mounted on a wire passing vertically through the 
camera at its center. The film is placed in a 
cylindrical envelope and fastened securely about 
the outer side. For the purpose of thermal insula- 
tion, the camera is packed in cotton in a wooden 
box lined on the inside with tin foil. There are 
three small openings in the box, one on the side 
for the slit system of the camera, one for the 
thermocouple and ground leads and the final one 
on top for a gas inlet. The outer end of the slit 
system is flush with the outside edge of the box 
to permit as great insulation as possible. 


TEMPERATURE MEASUREMENTS 


A copper-advance thermocouple was con- 
structed from No. 30 insulated copper wire and 
No. 30 insulated “‘advance”’ wire. Three advance 
wires and one copper wire lead to the variable 
junction, the wires at these points being soldered 
by means of tin solder. 

The potentials generated by this thermocouple 
were measured by means of a Leeds and Northrup 
Type K Potentiometer No. 169339 with lower- 
limit voltage readings at 0.000001 volt. The po- 
tentiometer was standardized by means of a 
Weston Standard Cell having a voltage of 1.0188. 
The same standard cell was used throughout all 
the experiments. 

Calibration of the thermocouple was carried 
out according to the method described by 
Adams." This method involves the determination 
of a correction curve showing deviations of e.m.f. 
from a set of standard values of e.m.f. as given by 
a general equation. The deviations are obtained 
by observing the potential generated when one 
junction is in ice water (0°C or 273°K) and the 
other junction is in a definite and accurate 
temperature bath. The standard values for the 

*Siegbahn, Spectroscopy of X-rays, Oxford University 
Press (1925). 

°N. W. Taylor, Rev. Sci. Inst. 2, 751-55 (1931). 

* Adams, J. Am. Chem. Soc. 36, 65 (1914). 
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low temperature calibration were obtained from 
Volume 1, page 57 of the International Critical 
Tables. The fixed temperature points used were 
(1) the boiling point of Oo, 


t= (—183.00+0.0126(p— 760) 
— 0.0000065 (p — 760)?) 


where p is the atmospheric pressure in milli- 
meters and ¢ is temperature in degrees centigrade; 
(2) the melting point of ethyl acetate (—83.6°C); 
(3) the melting point of mercury (—38.87°C) 
and (4) the melting point of carbon tetrachloride 
(—22.9°C).. The thermocouple was checked 
against a calibrated Bureau of Mines platinum 
resistance thermometer (Wheatstone Bridge, 
Leeds and Northrup No. 88131) and satisfactory 
agreement was obtained. The sensitivity of the 
apparatus used makes thermocouple readings ac- 
curate to +0.1°C. 


GENERAL PROCEDURE 


The general set-up is shown in Fig. 1. A and B 
are two 500 cc wide mouth bottles containing 
solid KOH for removal of COs, C is an 800 cc 
wide mouth bottle containing anhydrous calcium 
chloride for removal of HO, D is a 1/2 liter 
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Fic. 1. Low temperature x-ray camera. 


Dewar in which is placed a set of 1" 8” liquid 
air traps arranged in series, E is a 1-liter Dewar 
containing a meter of 6 mm copper tubing in the 
shape of a coil, F is a wooden box lined with tin 
foil. Inside this box, packed in cotton, is the 
camera G. K represents the thermocouple leads 
to the variable junction which is fastened to the 
wire specimen holder J, immediately below the 
point at which the x-rays coming through the 
slit H strike the specimen L. JV is a 250 cc Dewar 
which contains the ice and water mixture for cold 
junction. The specimen was ground to powder in 
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an agate mortar and placed on the wire specimen 
holder by means of collodion. Wood’s metal was 
used in connecting the thermocouple to the speci- 
men holder, since its melting point is so low that 
it can be applied easily and also will not affect 
the solder used in making copper-advance con- 
nections. A relatively large amount of Wood’s 
metal was used so as to make a thermal reservoir 
near the specimen. 

When the film is in position the only openings 
remaining in the camera are the tube on the 
upper side for gas inlet, the slit for x-ray entrance 
at the side and the tube at the bottom for en- 
trance of thermocouple leads. Constant tempera- 
ture was obtained by streaming cooled gas 
through the camera at such a rate as to counter- 
act exactly heat increase from without. 

The thermal capacity of the camera is rela- 
tively large so that small variations in the rate of 
flow of the gas will not be noticeable. However, 
during the period of an exposure, from two to 
three hours, air pressure and room temperature 
changes may require that slight adjustments be 
made in the rate of flow. Temperature deviations 
from the constant value desired are considered as 
being confined to a range of +0.5°C. 

Pictures of each of the compounds were taken 
at various temperatures throughout the range. 
Each is listed in a tabulated form below. In gen- 
eral, a picture of each substance was taken (1) at 
room temperature, (2) at the temperature cor- 
responding to the peak of the anomaly in the 
specific heat curve and (3) and (4) at tempera- 
tures above and below the anomaly. 


MATERIALS 


In these experiments we are dealing with the 
same material as used in the specific heat deter- 
minations. We are grateful to Dr. Charles G. 
Maier of the Bureau of Mines Pacific Experiment 
Station, Berkeley, California who has kindly sup- 
plied the compounds used. Complete details as to 
the preparation of these compounds are given in 
the reports on specific heat data. 


MEASUREMENTS 


The films were developed and measurements 
were made on negatives. The distances between 
outside edges of corresponding lines were meas- 
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ured to an accuracy of 0.1 mm. Under the condi- 
tions of the experiment it is impossible to give the 
exact dimensions of the camera. The tempera- 
ture measurements made are those of the central 
wire and not of the camera and hence the radius 
of the camera cannot be assumed. It is, there- 
fore, necessary to use some standard whose lat- 
tice constant is known so that the film may be 
calibrated. In relation to the compounds studied, 
the standard chosen must not only give lines 
which do not coincide with those of the compound 
but must also give enough lines throughout the 
range of the film so that calibration will be com- 
plete. No measurements of lattice constants of 
elements for low temperatures are recorded in the 
literature but, assuming that coefficient-of- 
expansion data are applicable to the side of a 
unit cell, we can calculate, on the basis of lattice 
constant at room temperature and coefficient of 
expansion data at low temperatures, the value 
of the lattice constant of the standard at the low 
temperatures. 

Mo was chosen as a standard for the com- 
pounds MnS, FeO and Fe;Q,. In the case of MnO, 
platinum wire was used, since the diffraction lines 
of MnO and Mo are coincident to the extent that 
they cannot be distinguished on the film. 

The coefficient of linear expansion of Mo wire is 
given by Schad and Hidnert™ as 


L,=Lo(1+5.15¢ X 10-§+0.00570# x 10), 


t being the temperature in degrees centigrade. 
This value agrees very well with results pre- 
sented by Disch." 

Coefficient-of-expansion data for platinum 
wire given by Dorsey were plotted and the 
average value over the range in question was 
obtained from the graph. 

Approximate values of the angle of reflection, 
26, and of sin? 6, were obtained by assuming 57.85 
mm as the effective diameter of the camera. 
Then by using the calculated lattice constants for 
various temperatures, the true values of sin? @ for 
the reflections from the different planes in the 
standard were calculated from Bragg’s equation. 


nd = 2d sin 8, 


12 Schad and Hidnert, Phys. Rev. 13, 148 (1919). 
13 J. Disch, Zeits. f. Phvsik 5, 173 (1921). 
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TABLE I. Data: MnO series: 111. MnO on Pt wire at 115°K. Film 266: Fe rad., 30 k.v. 15 m.a., 2 hr.: Lattice Constant 
Pt=3.904A, s=0.06 cm. 











Relative 2d—s Observed Calculated sin? 6 Corrected 

intensity (cm) 20 sin? 0 for Pt lines sin? 0 Indices sin? 0 
st. 12.34 122.21 0.7665 0.7674 (400) (16) (0.04795) 
m. 11.92 118.06 .7352 0.7362 (222) 
st. 11.10 109.93 .6704 .6748 (211) 
st. 9.89 97.95 5691 .5750 (222) (12) ( .04791) 
st. 9.32 92.31 .5209 .5275 (311) (11) ( .04795) 
m. 8.90 88.15 .4838 .4908 (220) 
v. Ss. 7.66 75.86 .3779 3836 (220) (8) ( .04795) 
m. 5.94 58.83 .2412 .2454 (200) 
st. 5.20 51.50 .1887 .1918 (200) (4) ( .04795) 
v. st. 4.47 44.37 .1419 .1439 (111) (3) ( .04796) 


sin? @=0.04795 (h?+k?+/*)---FeKa 
amno = (1.934)?/(4) (0.04795) =4.416A 








which, in case of a cubic lattice, may be written 
sin? @= (d?/4a*) (W?+2+P). 


The difference between the true and the ob- 
served values as plotted on a graph showing the 
relation between observed sin? @ and correction 
(difference between true and observed sin? 6). By 
making use of these correction curves, the correct 
values of sin* @ for all lines were obtained. Since 
the lines of the standard cover the range of the 
whole film, these correction curves also eliminate 
errors that may be due to the film itself, such as 
warping. The results are listed in tables showing 
the identification of lines with their respective 
crystal planes. Although lines due to Fe beta-rays 
were identified, only values resulting from lines 
due to Fe alpha-rays were used in calculating the 
lattice constants. The lattice constants are accu- 
rate to +0.002A. 

A typical set of data is presented in Table I, 
but on account of space similar data for the 
seventeen studied cases will-not be given. 

The lattice constants of the several substances 
are summarized in Table IT. 

The variations of the lattice constants with 


TABLE II, Lattice constants in Angstrom units at low 
temperatures. 











MnO °K Fe;O, °K MnS °K FeO °K 


4.436 299 8.363 299 5.210 299 4.290 299 

4.409 160 8.357 160 5.204 160 4.286 200 

4.416 115 8.363 114 5.204 143 4.284 186 

4.419 104 8.363 104 5.197 138 4.283 160 
5.192 130 














temperature are clearly anomalous. Note the 
minimum in the case of MnO and of Fe;O, at 
160°K. Refer to Figs. 2, 3, 4, and 5 for a plot of 
the lattice constant data. It should be stated that 
no change in crystal structure type appeared in 
the x-ray patterns. The only changes appear to 
be the irregularities in the expansion behavior. 


DISCUSSION 


Okamura" has carried out dilatometric meas- 
urements on magnetite. According to the graph 
given in his report Fe;0, expands with decrease 
in temperature from 173°K to 73°K. In the same 
report, it is stated that x-ray pictures by Nishi- 
yana revealed no change in structure. However, 
no description of the method used and no values 
for calculation are given. Two other physical 
phenomena of interest that occur at the transi- 
tion point are: (1) an abnormally great increase 
of specific resistance with decrease in tempera- 
ture and (2) a sudden fall in intensity of mag- 
netization after having gradually increased with 
decreasing temperature. 

Quoting from Okamura’s work: 


So far as the relative position of iron atoms in 
magnetite is concerned, alpha and beta magnetite have 
no noticeable difference. It is, however, to be noted 
that it is difficult to find the position of oxygen atoms 
in the iron lattice by Debye-Scherrer’s method of 
x-ray amalysis.... Thus the transformation of 
magnetite alphas@beta is probably due to a small 
displacement of oxygen atoms in the magnetite. 


144T,. Okamura, Science Reports Tohoku Imp. -Univ., 
First series 21, 2, 231-241 (1932). 
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Fic. 2. MnO: Relation between heat capacity, lattice con- 
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Fic. 3. Fe;04: Relation between heat capacity, lattice con- 


stant and temperature. 


Hence, in order to know the nature of this transfor- 
mation it is necessary to compare the Laue spots of 
alpha and beta magnetites. So far as any change of 
physical properties at the transformation temperature 
is concerned, magnetite shows a similar behavior to 
that of quartz. That is, quartz (SiOz) has an allotropic 
transformation at 573°; this change can be traced by 
those of thermal expansion, specific heat, electrical 
conductivity and magnetic susceptibility. Rinne found 


ELLEFSON AND N. W. TAYLOR 


7 MOL 


AL. 


¢ 


HEAT CAPACITY 


Fig. 


HEAT CAPACITY 


FIG. 


TEMPERATURE DEGREES KELVIN 


4. MnS: Relation between heat capacity, lattice con- 
stant and temperature. 
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5. FeO: Relation between heat capacity, lattice con- 
stant and temperature. 


by means of x-ray analysis that the transformation of 
quartz alphazbeta is due to a small displacement of 
oxygen atoms. Thus in all probability the trans- 
formation in magnetite is similar to that in quartz. It 
is very interesting to note that the transformation in 
magnetite and in quartz which consists only of the 
small change of distribution of the oxygen atoms in 
the iron lattice affects the electrical and magnetical 
properties to a considerable degree. This marked 
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effect is probably caused by some change in the 
valence electron in the atoms, which may reasonably 
take place as the result of small displacement of 
atoms in that lattice. Thus the transformation of a 
substance may consist of three changes, that is, 


(1) change of lattice, 
(2) change of state of electrons in atoms, 
(3) change of orientation of atoms, 


In the case of magnetite and quartz, change (1) is 
probably small but (2) and (3) are considerable. 


By dilatometric measurements carried out by 
Simon and Bergmann!’ it is shown that am- 
monium chloride increases in volume at the 
transition point, while ammonium bromide and 
ammonium dihydrogen phosphate gradually 
contract during the transition interval. Kracek'® 
noted an abnormal expansion of sodium nitrate. 

Investigations by Southard, Milner and Hend- 
ricks!? on normal amyl ammonium chloride indi- 
cate a similar behavior. 

On the basis of a wave-mechanical treatment, 
Pauling'® shows how many gradual transitions 
from an oscillatory state to one of rotation can 
arise in solids where the moments of inertia are 
small and the intermolecular forces are also small. 
Such transitions usually are ‘“‘foreshadowed on 
the low temperature side by an abnormal in- 
crease in heat capacity.” 

S. B. Hendricks, E. Posnjak and F. C. Kracek 
have investigated ammonium nitrate’ and ex- 
plain some of their results on the basis of rotation 
of nitrate ions. The abnormalities in ammonium 
halides to which reference was made at the begin- 
ning of this paper are explained on the basis of 
rotation of ammonium ion. 

On the other hand, in solids having strong in- 
termolecular ionic or atomic forces such transi- 
tions from oscillation to rotation are less likely. 
In view of the fact that MnO, FeO, MnS and 
Fe;O, are all rather polar compounds of high 
melting point, hardness, etc., they would seem to 
fall into this latter class. We must therefore look 


%F, Simon and R. Bergmann, Zeits. f. physik. Chemie 
B8, 255 (1930). 

* F.C. Kracek, J. Am. Chem. Soc. 53, 2609-24 (1931). 

Southard, Milner and S. B. Hendricks, J. Chem. 
Phys. 1, 95-162 (1933). 

*L. Pauling, Phys. Rev. 36, 430 (1930). 

”S. B. Hendricks, E. Posnjak and F. C. Kracek, J. Am. 
Chem. Soc. 54, 2766-86 (1932). 


for another interpretation of their anomalous 
specific heat and expansion behavior. 

As far as behavior of the specific heat curve is 
concerned, Fe at 765°C is quite analogous to the 
compounds we have studied. Since this ab- 
normality was explained by magnetic changes we 
may expect them in our compounds. Okamura’s 
work bears out this argument. Of interest in this 
connection is the fact that CuO shows abnormal- 
ities in its specific heat behavior while CuO does 
not.! Cu,O is diamagnetic whereas CuO is para- 
magnetic. All of the compounds which we have 
studied are paramagnetic. J. H. Van Vleck” 
states: 


Even in paramagnetic bodies some anomalies in 
specific heat curves should be expected. There is a 
considerable amount of spin coupling between 
paramagnetic ions due to the Heisenberg exchange 
effect, provided the material is one of fairly high 
magnetic concentration, e.g., MnS, MnO, Fe;QO.. 
This is presumably the cause of anomalies in MnS and 
CuO. It might be described as a potentially ferro- 
magnetic body above the Curie point. It should alter 
the specific heat but it is not at all clear that it should 
give a cusp in the curve. In fact, the conventional 
theory gives an anomaly only at the Curie point with 
no change in specific heat above. 


The specific heat curves of MnO and FeO; 
show striking similarity. The abnormality in each 
case is very great and is confined to about 20° 
and the peak is at practically the same tempera- 
ture (115°K for MnO and 114°K for Fe;Q,). 
Above the abnormal range and up to 300°K, the 
C, curve is a continuation of the initial stage and 
also is concave to the temperature axis in each 
case. Curiously enough both substances show 
contraction as the temperature rises through the 
abnormal range from 104° to 160°, followed by 
expansion at higher temperatures. It may be of 
interest to bear in mind that Fet** ion in mag- 
netite has essentially the same electron structure 
as Mn+ in MnO. In the case of MnS, the bond 
is probably not of a simple ionic type because of 
polarization of the sulfide ion. 

The specific heat curves of MnS and FeO are 
similar to one another in shape but different from 
MnO and Fe;0,. They both have relatively large 
expansion through the abnormal interval fol- 


20 J. H. Van Vleck, private communication. 


. 








64 B. S. ELLEFSON AND N. W. TAYLOR 


lowed by smaller expansions at higher tempera- 
tures. 
This is seen from the following data: 


Temperature range 160—200°K 200-299°K 
Coefficient of expansion of FeO 1.75X10~° 0.94 x 10-5 
Temperature range 130-160°K 160-299°K 


Coefficient of expansionof MnS_ 7.7 X10~* 0.83 x 10-5 


It is difficult to say that there is a real differ- 
ence in the behavior of MnS and FeO as com- 
pared with MnO and Fe;O, but the data would 
seem to indicate it. Magnetic susceptibility 
measurements over the abnormal temperature 
range of each of the substances may bring out 
more clearly any difference. 
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The Calculation of Thermodynamic Quantities from Spectroscopic Data for Poly- 
atomic Molecules; the Free Energy, Entropy and Heat Capacity of Steam 


A. R. Gorpon, Chemistry Department, University of Toronto 
(Received December 11, 1933) 


The methods already developed by Gordon and Barnes 
for computing thermodynamic quantities from the band 
spectrum of a diatomic molecule are extended to deal with 
the case of the triatomic unsymmetrical rotator. The free 
energy, entropy and heat capacity of steam are computed 
for the temperature range 298.1-1500°K from the spectro- 
scopic data of Mecke and his associates. From the free 


energy so calculated, the known spectroscopic free energies 
of carbon monoxide, carbon dioxide and hydrogen, and the 
thermal data of Rossini, revised values of the equilibrium 
constant for the water-gas reaction are computed. The 
numbers so obtained are in close agreement with the 
experimental results of Neumann and Koehler and of 
Emmett and Schultz. 





HE calculation of thermodynamic quan- 
tities from the band spectrum of a diatomic 
molecule has been discussed in detail by various 
authors! and a number of procedures have been 
developed which avoid the principal difficulties 
of summation. The corresponding calculation for 
the polyatomic case has, in general, received 
scant attention primarily, no doubt, because the 
requisite experimental data have been lacking. 
In a recent series of papers, however, Mecke and 
his associates? have presented a reasonably satis- 
factory analysis of the rotational-vibrational 
spectrum of steam and although there are still 
some uncertainties as to the changes in the spec- 
troscopic constants in the higher vibrational 
states, it is now possible to obtain reliable values 


of the free energy, entropy and heat capacity for 
moderate temperatures at least. Moreover, the 
series which arise in this case are typical of those 
which will occur in similar calculations for other 
triatomic unsymmetrical rotators, e.g., SO. and 
H,S, and it therefore seems justifiable to consider 
the case of steam in some detail. Indeed, steam 
provides a particularly good test for any method 
of approximation, since its rotational-vibrational 
interaction terms are exceptionally large. 

The approximations discussed here are an 
extension of those developed by Gordon and 
Barnes’ for diatomic molecules and the problem 
is really one of reducing the ‘‘state sum”’ and its 
derivatives to such a form that the tabulated 
diatomic functions can be used in their evalu- 
ation. 





THE EVALUATION OF THE STATE SUM AND THE CALCULATION OF THE FREE ENERGY 


The rotational-vibrational energy levels for a triatomic unsymmetrical rotator will, in general, 
depend on the three fundamental frequencies of the molecule and on its three principal moments of 
inertia. The frequencies are distinguished by the subscripts 6, z and o, the corresponding vibrational 
quantum numbers being v5, v, and v,; in the expression for the energy levels, not only will terms 
occur involving the first and second powers of the vibrational quantum numbers but also ‘‘cross 
terms’’ of the type vs‘v,, etc.; moreover, the moments of inertia A, B and C will themselves be 





1 Mulholland, Proc. Camb. Phil. Soc. 24, 280 (1928); 
Sutherland, Proc. Camb. Phil. Soc. 26, 402 (1930); 
Giauque and Overstreet, J. Am. Chem. Soc. 54, 1731 
(1932); Gordon and Barnes, J. Chem. Phys. 1, 297 (1933); 
ay Phys. Rev. 43, 364 (1933); J. Chem. Phys. 1, 576 
1933), 


* Mecke, Zeits. f. Physik 81, 313 (1933); Baumann and 
Mecke, Zeits. f. Physik 81, 445 (1933); Freudenberg and 
Mecke, Zeits. f. Physik 81, 465 (1933). 

* Gordon and Barnes, J. Chem. Phys. 1, 297 (1933). 

* Throughout this paper it is assumed that A <B<C. 
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functions of v3, v, and v,. If the molecule is of the type RX Y with non-spinning nuclei, the sum- 
mation over the rotational levels for a given vibrational state is replaced by its limiting value for 
high temperatures: 







Qr~exp (o8c/4)-Qv- (1 +oR8¢/12—opo¢/1204), (1) 





where 


qv =82°M ,(2rkT/h?)' = (2/cacgoc)}; 












M,=(A,B,C,)}; 
oA, OB, Oc =(h?/8x*RT)-(1/A>d, 1/Bo, 1/C.); opc=(opec)! 








and the subscript v stands for v3v,v,. The error introduced by this approximation is discussed in the 
note at the end of this paper; it is there shown that the approximation is justifiable in any case likely 
to occur in practice, provided the temperature is 300°K or more. 

If, as in the case of steam, there are two equivalent atoms in the molecule each possessing js 
units of nuclear spin, the rotational state sum becomes 


Qrexp (opc/4) ‘dv’ (1+¢p¢/12 —opoc/12c4) . (2ist+ 1)?/2. (2) 








The vibrational terms in the expression for the energy levels will be of the type 


ey = (€vst (Ev) et (€vdotei, 










(3) 







where 






(€v)s = (Vs t+ 2)as— (vs +3)? wsx;, 


€i = — Vor VV er — Varo UVa — Vad Vos 






Thus the rotational-vibrational state sum, Qry, 







with definitions similar to (4) for (€,), and (€y).. 
may be written 


Qrv =z exp [—e,/kT ]-exp (oBc/4)-qvX At+oenc/12—opo¢/1204) a (2js+1)*/2, (6) 










where the summation is to run over all values of v5, v, and v.. 

Since A,, B, and C, can all be expressed in terms of their values in the lowest vibrational state 
plus certain linear and quadratic terms in the vibrational quantum numbers, it follows that q, can 
be replaced by an expression of the type 


v=o 1 t+asstaWstaerotbwet:--+C5nViet"**), (7) 













where gp is the value of g, for v5 =v,=v,=0, the a’s, b’s and c’s are numerical constants and the dots 
after a plus sign indicate two terms obtained by permuting 6, 7 and o in the term immediately 
preceding it. 

Since ¢;/kT is in general small, exp [e;/kT_] can be expanded in a power series and by an obvious 
rearrangement the logarithm of the state sum can then be written 


In Qrev+Eo/RT = (In Qvy+e0/kT)s+ (In Qo +¢0/kT) x +(In Qu t+e0/kT) 
+opc/44+In (2j7s+1)%qo/2+In (1+7;), (8) 
















where 
r= (yor/RT) + (5) 5(d) e+ ++ +3 (Yar/kT)*- (o> (P) ate ++ +4(yex/kT) + (yee/kT)-(B)a(0*) «(Bo 
f+++fag(d)st- ++ +b3(v%)st+--++c52(0)3- (5) +°**+0Bc/12—ep0c/1204, (9) 
(In Qy+eo/kT)s=In 2 exp [—(vsws— (vs? +03) wsxs)/RT J, (10) 
(5)s=(Zvs-exp [— (ev)s/kT ])/(Qv)s, (11) 
(v?) 5 = (Zvs?-exp [— (€o)s/kT])/(Qv)s, (12) 
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with similar definitions for symbols involving + and o. Since numerical values of (In Qy+eo/kT), 
5 and v2 can be found from Gordon and Barnes’ Tables I, II, V and IX or by means of the approx- 
imations discussed in their Example 4, for given values of w= w/kT and x, it follows that the numer- 
ical value of the expression on the right of Eq. (8) can be found at once. If 7; is small in comparison 
with unity, Eq. (8) can be still further simplified, giving 


In Qev+Eo/RT = (In Qy+e0/kT)s+: ++ +1n (275+1)*Go0/2 
+ (ysr/RT)+ (ds: (Bat+ ++ +as(B)st+ ++ +0s(0*)st+onc/3—exoc/1204. (13) 


THE CALCULATION OF THE ROTATIONAL-VIBRATIONAL ENTROPY AND HEAT CAPACITY 
Since Sev=R In Qrev+ RTO In Qev/dT and Cry =TdSry/dT it follows from Eq. (8) that 
Srv/R=(Sv/R)s+(Sv/Rv)s+(Sv/R)etln (2js+1)%o/2+3/2+In (1+1r45) +(Tari/dT)/(1+r:), (14) 
Crev/R=(Cy/R)st+(Cv/R)r+(Cv/R)o+3/2 
+ {Ta(Tdr,/dT)/dT} /(A+r:) +(Todr:/dT)/(A+r:) —(Ter,/dT)?/(1+r,)*, (15) 


where (Sy)s and (Cy), stand for the entropy and heat capacity of an anharmonic oscillator whose 
energy levels are given by Eq. (4), with similar definitions for symbols involving 7 and o. The 
expressions for 7'r;/8T and Td(Tdr;/dT)/dT can be readily written down by means of the typical 
relations, 


T06/dT =s,—3, (16) 

Tdv/dT =s.—0, (17) 

Tas,/dT =e, (18) 

T s2/AT =Cs, (19) 

T(0) 5: (5) x/OT = (S1)s° (0) e+(B)s- (51) —2(5)5° (8), (20) 
T9{(d)s- (0) x/T} OT = {(s1)s- (0) e+ (05° (S1)¢—3(0)5- (0) x} /T, (21) 
T0(S1)a- (0) x/OT = (61) 5° (0) 2 +(S1) 8° (S1) x — (Si) 8° (0) x, (22) 
TO{(si)s- (6) 2/T}/OT = { (C1) 3(0) «+ (S1)8- (51) e—2(S1) 8: (0) x} /T. (23) 


Since numerical values of Sy/R, Cy/R, 51, S2, ¢: and ¢z can be found from Gordon and Barnes’ 
tables for given values of w/kT and x, Eqs. (14) and (15) may be used to find an approximate value 
of the rotational-vibrational entropy and heat capacity. Alternatively and to a somewhat poorer 
approximation, from Eq. (13) 


Srv/R=(Sy/R)s+(Sy/R) «+(Sv/R)o+ln (2j7s+1)?qo/2+3/2 
+ (ysx/RT) + {(S1)8- (0) +(O)s- (S1)e—2(d)s-(H) e}+--- 
+45($1)s+@+(S1)e¢+Ge(S1)o+bs(se)st+:-- (24) 


Cev/R=(Cy/R)st+(Cy/R)2+(Cy/R)e+3/2 
+ (yor/RT) + { (1) 5: (8) e+ (5° (C1) x +2(51) 5° (51) —4(51) 5° (8) x —4(0) 8- ($1) x + 6(8) 5° (0) x} 
+ +++ +05(C1)stax(C1)e+e(C1)e+bs(C2)st+-+-. (25) 
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TABLE I. T=1500°K, he/k=1.43242. 































6 T o 6 T o 
w=w/kT 1.5615 3.5715 3.6613 v 0.287 0.033 0.028 
x 0.01223 0.01872 0.01017 Si 0.859 0.153 0.134 
In Qy +60/kT 0.2356 0.0285 0.0260 1 1.398 0.446 0.403 
+0.0136 +0.0044 +0.0021 a 
Sy/R 0.6504 0.1318 0.1226 v 0.46 —_ — 
+0.0327 +0.0161 +0.0080 
Cy/R 0.8194 0.3797 0.3629 S2 1.73 — _- 
+0.0374 +0.0317 +0.0164 
C2 4.6 — — 




















THE FREE ENERGY, ENTROPY AND HEAT 
CAPACITY OF STEAM FOR TEMPERATURES 
uP To 1500°K 


From the data of Freudenberg and Mecke,? 
the vibrational constants in cm™ are: w;=1635, 
@®,=3740, we= 3834, wsxs=20, wexr=710, Wexe 
= 39, ys,= 20, yxe=109, y.5=20. The moments 
of inertia are 


10° X A =0.996+0.045v,+0.026v, —0.098v;, 
10” x B= 1.908 + 0.014v, +0.033v, —0.034v5, (26) 
10° X C=2.981+0.047v,+0.062v,+0.06205. 


While there may be some doubt as to the accu- 
racy of Freudenberg and Mecke’s expressions for 
the moments of inertia in the higher vibrational 
states, their equations seem reasonably accurate 
for the states which are of importance in cal- 
culations for temperatures below 1500°K. From 
Eqs. (26) 


M,=2.3801 X 10-©(1+0.0342v,+0.03210, 


— 0.04770; —0.00147»,?), (27) 


the other quadratic terms in the expression for 
M, making negligible contributions to the cal- 
culated thermodynamic quantities for the tem- 
perature range considered. 

Table I gives (for 7=1500°K) the values’ of 
w/kT and x for the three fundamental fre- 
quencies and also the values of the functions 
necessary for the calculation; these latter were 
obtained by interpolation in Gordon and Barnes’ 
tables. 


5 The universal constants used in computing Tables I and 
II were those given in International Critical Tables, Vol. I: 
h=6.554X1077,  k=1.372K10-", C=2.998610" R 
= 1.9869; hence the additive constant for the translational 
free energy equation is —7.2671. 











For this temperature, In go/2 = 6.19597, ogc/3 
=().00369, ogo¢/1204=0.00039. Hence, exclud- 
ing the nuclear spin contributions to the free 
energy and entropy, the approximate values of 
In Qevt+Eo/RT, Srv/R and Crv/R from Eggs. 
(13), (24) and (25) are 6.4975, 8.6261 and 3.113, 
respectively; the “‘exact’’ values® from Eq. (6) 
are 6.9474, 8.6258 and 3.112. 

The contributions of the vibrational inter- 
action terms involving ys, Yro and yes to 
In Qev+Eo/RT are only 0.00018, 0.00009 and 
0.00015 for this temperature; the corresponding 
items for Sry/R are 0.00102, 0.00071 and 
0.00089, while for Cry /R they are 0.0039, 0.0039 
and 0.0035, respectively. The y-terms are thus 
of much less importance than the rotational- 
vibrational interaction terms, i.e., those involving 
a3, Gx, a@¢ and 63, which contribute —0.01236, 
—0.03403 and —0.0453 to In Qrev+Eo/RT, 
Srv/R and Crv/R, respectively. 

Table II gives the numerical values’ of 
—(F°—E,°)/T, S° and Cp° for steam for the 
temperature range 298.1-1500°K, the rotational- 
vibrational contributions being computed by 
means of Eqs. (13), (24) and (25). None of the 
values of —(F°—£,°)/T and S° includes the 
nuclear spin contribution of the two hydrogen 
atoms, 2R In 2=2.754 cal./deg.; this must be 
added to the entries if the table is to be used in 
conjunction with Giauque’s table for hydrogen.’ 
The third decimal place has been retained in the 
table merely for convenience in interpolation. 

From Table II and the data for hydrogen’ and 































6 These values are exact in the sense that the only error 
is that introduced by the use of Eq. (1); in the note at the 
end of this paper it is shown that this error is negligible for 
1500°K. 

7 Giauque, J. Am. Chem. Soc. 52, 4816 (1930). 
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TABLE II. The free energy, entropy and heat capacity of steam. 








— (F°— Ey°)/T S° C°p 


—(F°—E,°)/T a 





8.000 
8.002 
8.066 
8.155 
8.260 
8.379 
8.504 
8.635 
8.771 
8.910 
9.053 
9.199 


37.179 
37.230 
38.452 
39.513 
40.452 
41.296 
42.062 
42.765 
43.415 
44.020 
44.587 
45.121 


45.101 
45.151 
46.389 
47.472 
48.439 
49.315 
50.119 
50.864 
51.561 
52.216 
52.836 
53.425 





53.987 
54.525 
55.043 
55.542 
56.023 
56.489 
56.941 
57.380 
57.807 
58.223 
59.022 
59.783 


45.627 
46.106 
46.563 
46.999 
47.418 
47.820 
48.206 
48.579 
48.940 
49,289 
49.956 
50.586 








oxygen,® AS°oo.1 for the reaction H2+O2/2 
=H,0(g) is — 10.64 cal./deg.; from this in com- 
bination with Rossini’s accurate value? of AH95.1, 
viz., —57,809 cal., AF°oos.; for steam is — 54,637 
cal. For liquid water, S°2o9s.1 is 16.75 cal./deg. 
and AF° 298.1 is — 56,690 cal. 

The data of Table II can also be combined 
with the spectroscopic free energy data for 
hydrogen,’? carbon monoxide” and carbon diox- 
ide" and the thermal data of Rossini® ” to 
compute the equilibrium constant for the water- 
gas reaction CO+H2O=CO2+He. AS°298.1 is 
here —10.09, AHogs.; is —9814; hence AF°a9.1 
is —6806 cal. and AE,° is —9646 cal. Table III 
gives the resulting values of RlnK=RIn 
(Pco' Pu,o/Pco,*Pu,) for the range 600-1500°. 
The numbers in Table III are on the average 
about 0.07 cal./deg. less than those previously 
computed" and for the range 900—-1259° agree 
even more closely with the results of Neumann 
and Koehler’s experiments;' for example, Neu- 


TABLE III. The equilibrium constant for the water-gas 
reaction. 








1000 
—0.69 
1500 
+1.88 


900 
—1.62 
1400 
+1.53 


T°K 800 
—2.81 
1300 


+1.12 


600 

— 6.57 
1100 
+0.04 


700 
—4.40 
1200 


Rink +0.63 








‘Johnston and Walker, J. Am. Chem. Soc. 55, 172 
(1933), 

* Rossini, Bur. Standards J. Research 6, 1 (1931). 

Clayton and Giauque, J. Am. Chem. Soc. 54, 2610 
(1932). 

"Gordon, J. Chem. Phys. 1, 308 (1933). 

® Rossini, Bur. Standards J. Research 6, 37 (1931). 

3 Neumann and Koehler, Zeits. f. Elektrochemie 34, 218 
(1928), 


mann and Koehler found K =1.61 for 1259°K, 
while from Table III, R In K for this temperature 
is 0.93, i.e., K=1.60. Their results for 15 tem- 
peratures in this range show a mean deviation 
from the numbers obtained by interpolation in 
Table III of less than +0.02 cal./deg. in R In K; 
for the range 700°-900°, however, the mean 
deviation is somewhat larger—about +0.14 in 
Rin K. The indirect measurements of this 
equilibrium by Emmett and Schultz" also agrees 
satisfactorily with Table III; with three excep- 
tions (773°, 823°, 923°) their results for 8 tem- 
peratures between 723° and 1073° agree with 
the data of the table within +0.04 cal./deg. in 
Rin K. 


EVALUATION OF THE ROTATIONAL STATE SUM 
FOR A RIGID UNSYMMETRICAL ROTATOR 


The procedure that has been adopted almost 
universally in the past when it has been necessary 
to evaluate the rotational state sum for a rotator 
with three different moments of inertia has been 
to replace the sum by its limiting value for high 
temperatures 


(28) 


Or=(14/caonoc)', 


where oa, og and o¢ are as defined under Eq. (1) 
and suitable symmetry and nuclear spin factors 
are to be applied when necessary. While Eq. (28) 
is a valid approximation at very high tem- 
peratures, there is no general rule for estimating 
its error in any actual case. In the special case 
that two of the three moments of inertia are 


4 Emmett and Schultz, J. Am. Chem. Soc. 51, 3249 
(1929); 52, 892 (1930); 55, 1376 (1933). 








equal, Viney" and Kassel! have deduced a rela- 
tion which serves to evaluate Qe to a higher 
degree of accuracy than the uncorrected Eq. 
(28). If the moments of inertia are A, C and C, 
their equation may be written, to a first order of 
approximation 


Qr=exp (¢¢/4)-(/o40°%c)! 
-{1+0¢/12—07¢/12044+---}. (29) 


In the general case, however, where A +B+C, 
their method cannot be applied, since the energy 
is not an explicit function of the moments of 
inertia. It is apparent, however, that the ex- 
pression on the right of (29) for a rotator with 
moments of inertia A, (BC)! and (BC)!, should 
be a better approximation to the value of the 
actual rotational state sum than the uncorrected 
expression on the right of (28). On this plausible 
assumption the rotational state sum is replaced 
by the expression on the right of Eq. (1). 

To obtain an estimate of the error actually 
introduced by the use of Eq. (1), the energy 
levels were calculated for the special case 
A=1xX10-, B=2x10-", C=3X10™, and the 
state sum and its derivatives were computed. 
These constants are roughly those for steam and, 
since they correspond to a molecule as unsym- 
metrical and with as small moments of inertia 
as is likely to be met with in practice, they 
provide a particularly good test of Eq. (1). 

The algebraic equations from whose roots the 
eigenvalues for the rotator can be obtained, have 
been tabulated by Dennison" and by Nielsen,” 
For values of J, the rotational quantum number, 
<8, the energy levels were found by solving 
their equations (see Table IV); for J$9, the 
asymptotic formulae of Kramers and Ittman’® 


% Viney, Proc. Camb. Phil. Soc. 29, 142 (1933). 
16 Dennison, Rev. Mod. Phys. 3, 280 (1931). There are 
two misprints in Dennison’s cubic for J=7; he gives 


x8 — 56ax*?+ (784a?+ 7768)x —2304a*—10584a8=0. (A) 


By a change of variable, Nielsen’s!’ corresponding cubic 
can be written in Dennison’s notation 


x3 — 56ax*+ (7840? + 7848)x —2304a*—11584a8=0. (B) 


Eq. (B) fulfills the condition that the energy levels should 
be the same whether a and £ are replaced by aif, a2 or 
383. 

17 Nielsen, Phys. Rev. 38, 1432 (1931). 

18 Kramers and Ittman, Zeits. f. Physik 58, 217 (1929). 
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were used. While the latter are not particularly 
accurate as far as the calculation of individual 
levels is concerned, it was found by actual com- 
putation for J=7 and J=8 that they give the 
contributions to Qr and its derivatives with an 
accuracy of about 1/5 of a percent. Since even 
for T =330°K (the highest temperature for which 
the calculation was made), the levels with J>8 
contribute less than 3 percent of Qe, 10 percent 
of Qr’ (see Eq. (33)) and 22 percent of Qp”’ 
(Eq. (34)), the error introduced by using the 
asymptotic formulae is, at most, a unit or two 
in the fourth decimal place of In Qe and Sr/R 
and in the third decimal place of Cr/R. For 
lower temperatures for which the relative con- 
tribution of the high J states is less, the uncer- 
tainty will be still smaller. 

For a rigid rotator with three different mo- 
ments of inertia, there are 2/+1 distinct levels,’ 
each (2J+1)-fold degenerate, for each value of 
the rotational quantum number J; these 2J/+1 
levels are distinguished by a quantity 7+ which 
takes the 2/+1 values —J, J+1, --- J—1, J, 
the lowest level (for the given J) being assigned 
7=-—J the next lowest r= —J+1 and so on. 
The rotational state sum thus becomes 


@ 


+J 
Or=dX DY (2I4+1)-e-esr!/*? (30) 


J=0 r=—J 


and the rotational entropy and heat capacity are 
given by 


TABLE IV. Values of (8x*e7,/h?) X10. 
A=1X10-", B=2x10-*, C=3X10-*°. 








J=0: 0. 

J=1: 5/6, 4/3, 3/2. 

J=2: 2.24267, 2.83333, 3.33333, 4.83333, 5.09067. 

J=3: 4.83333, 5.06325, 6.05848, 7.33333, 7.50000, 
10.27008, 10.27485. 

J=4: 7.88133, 8.00491, 9.63908, 10.63900, 11.09047, 


13.66176, 13.69425, 17.69433, 17.69487. 


: 11.5824, 11.6417, 14.0188, 14.7278, 15.6406, 17.9080, 
18.0308, 21.9389, 21.9436, 27.1170, 27.1172. 


J 
J=6: 15.934, 15.961, 19.126, 19.574, 21.118, 22.999, 
J 


ll 
w 


23.333, 27.053, 27.075, 32.207, 32.207, 38.540, 38.540. 


=7: 20.945, 20.955, 24.900, 25.153, 27.460, 28.916, 29.635, 
33.041, 33.122, 38.166, 38.170, 44.473, 44.473, 51.962, 
51.962. 


J=8: 26.616, 26.620, 31.310, 31.441, 34.603, 35.634, 36.930, 
39.899, 40.120, 45.003, 45.020, 51.276, 51.276, 58.743, 
58.743, 67.385, 67.385. 
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Sr/R=OQr'/Qrtln Qe, (31) 
Cr/R=Qr"’/Qr—(Qr’/Qr)’, 
Or’ = LU(2I+1)-(ex-/kT) exp [—ey,/RT ], (33) 


Qn! =22(2I+1)- (ese/kT)? exp [—es-/kT]. (34) 


(32) 


Table IV gives the numerical values of Qe, Q’r 
and Qr” as computed from Eqs. (30), (33) and 
(34) together with the resulting In Qg+Eo/RT, 
Sr/R and Cr/R for T =330.44°, 198.26°, 132.17° 
and 99.13°K, i.e., for 10-h?/8?kT =0.12, 0.2, 
0.3 and 0.4; the table also gives the values ob- 
tained by using the approximate expression (1) 
instead of the exact expressions (30), (33) and 
(34). It should be noted that if the uncorrected 
Eq. (28) had been used in place of Eq. (30), the 
resulting value of In Qz would have been 4.6486 
for 330.46° and 2.8427 for 99.13°. It is evident 
from the table that with the assumed moments 
of inertia and for room temperature and up, the 
use of Eq. (1) will not introduce an error of more 
than 0.001 cal./deg. in F°/T or S°, or of more 
than 0.002 cal./deg. in Cp°. Since the case chosen 
for illustration is a particularly unfavorable one, 
it seems safe to conclude that Eq. (1) may be 
used with confidence in similar calculations for 
other unsymmetrical rotators. For low tem- 
peratures, however, the calculated thermody- 
namic quantities deviate widely from those 
predicted by Eq. (1); for example, curve a of 
Fig. 1 shows the low temperature rotational 
heat capacity curve for the rotator of Table IV. 
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Fic. 1. Cr/R as a function of temperature for an un- 
a ac rotator A=1X10™, B=2X10™"%, C=3 
40 


In the previous discussion, all questions of 
symmetry and of nuclear spin have been ignored. 
If the rotator possesses two equivalent atoms 
with non-spinning nuclei and if, as in the case of 
steam, the symmetry axis of the molecule coin- 
cides with the intermediate moment of inertia, 
Dennison" has shown that the rotational levels 
for a given J are alternately symmetric and anti- 
symmetric in the rotational part of the wave 
function; for J/=0, 2, 4, ---, the lowest level is 
symmetric and for J=1, 3, ---, antisymmetric, 
i.e., for either odd or even J, the levels with 
7=0, +2, +4, - are symmetric and those 
7=+1, +3, ---, antisymmetric. In carrying out 
the summations of Eqs. (30), (33) and (34) only 
those rotational levels are to be counted which 
are of the required symmetry character and, for 
any temperature so high that the summations 
over the symmetric and antisymmetric states 
are the same, this is obviously equivalent to 
adding the usual Ehrenfest ‘‘symmetry term”’ 
—RIin2 to —F°/T and to S®° and leaving Cp° 
unaltered. For the rotator of Table IV at the 
highest temperature of Table V, this condition is 
fulfilled, e.g., the sums over the symmetric and 
antisymmetric levels in Qr are 52.963 and 
52.957, respectively, but for low temperatures 
this is not the case. For example, if the electronic 
and vibrational parts of the wave function are 
such that only symmetric rotational states are 
to be counted, Cr/R for 198.26°, 132.17° and 
99.13° is 1.508, 1.516 and 1.513, respectively; 
for still lower temperatures, the Cpr/R-tem- 
perature curve for such a symmetrical molecule 
is shown in curve 6 of Fig. 1. 

If the equivalent atoms, as in steam, each 
possess one-half unit of nuclear spin (the ex- 


TABLE V. A=1X10-*, B=2X10-", C=3X10-*, 








T°K 330.44 198.26 132.17 99.13 
Or, Ea. (30) 105.92 49.700 27.4189 18.0677 
Or’, Eq. (33) 157.34 73.290 39.9757 26.0384 
Or”, Eq. (34) 392.76 182.894 99.6460 64.7483 
In Qr+Eo/RT 

Eq. (30) 4.6627 3.9060 3.3112 2.8941 
In Or+Eo/RT, 

Eq. (1 4.6633 3.9068 3.3109 2.8916 
Sr/R, Eq. (30), 

33 6.1482 5.3807 4.7692 4.3353 
Sr/R, Eq. (1) 6.1486 5.3824 4.7742 4.3427 
Cr/R, Eq. (30), 

33), (34) 1.501 1.505 1.509 1.507 
Cr/R, Eq. (1) 1.500 1.500 1.500 1.500 
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tension to other js values is obvious) and if the 
symmetry requirement of the molecule as a 
whole is such that states symmetric in the rota- 
tion are to be combined with states antisym- 
metric in the spin and vice versa, the weight 
factors will be 2/+1 for the “‘para’”’ levels, i.e., 
those with r=0, +2, ---, and 3(2/+1) for the 
“ortho” levels, those with 7= +1, +3, ---. For 
high temperatures this is equivalent to multi- 
plying Qr, Qr’ and Qr’’, computed in the 
absence of symmetry and spin, by (1/2+3/2), 
that is, to adding the usual symmetry and spin 
terms —R1In2+2RIin2 to —F°/T and Sy, cf. 
Eq. (2). This method of allowing for symmetry 
and spin is thus justified for the rotator of 
Table IV at the highest temperature of Table V 
but is not valid for any temperature so low that 
the high temperature 3:1 ratio of ortho to 
para molecules does not hold (see Table VI 


TABLE VI. Percentage para molecules in the equilibrium 
mixture. 








T°K 330.44 132.17 66.09 39.65 26.44 
Percent 25.00 25.12 29.91 26.52 30.04 
T°K 15.86 11.33 7.93 5.66 4.96 
Percent 45.58 65.82 87.42 97.41 98.86 








which gives the equilibrium percentage of para 
molecules at various temperatures for the rotator 
of Table IV). The rotational heat capacity curve 
of such an equilibrium mixture is shown in c of 
Fig. 1; for higher temperatures, Cr/R approaches 
the limiting value 3/2 as in the case of curves 
a and bd; for example, Cr/R=1.503, 1.505, 1.503 
and 1.501 for 99.13°, 132.17°, 198.26° and 
330.44°, respectively. 

Since the ortho and para states are non-com- 
bining, it is probable that a mixture cooled from 
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high temperatures would retain at low tem- 
perature the 3 : 1 ortho-para ratio in a state of 
“frozen” equilibrium. The rotational heat capac- 
ity curve for such a metastable mixture is shown 
in curve d of Fig. 1; for T>50°K in the case 
chosen for illustration, curves c and d become 
practically identical. 

The possibility of a frozen equilibrium of this 
kind is of interest in connection with the ‘“Third 
Law”’ entropy of steam, which is nearly 0.9 
cal./deg. less than the spectroscopic value. If ice 
consists of such a 3 : 1 mixture and it is assumed 
that the rotational states of steam carry over 
into the crystal, the ortho molecules will be 
trapped at low temperatures in the lowest state 
antisymmetric in the rotation i.e., J=1, r= —1; 
this state has a multiplicity of 9 (3 from rotation 
and 3 from spin) and as a consequence the 
“zero-point” entropy over and above the high 
temperature spin entropy 2RIn2 will be 
3R/4-In 3=1.64 cal./deg. just as in the case of 
hydrogen. Giauque and Ashley” have pointed 
out, however, that at temperatures as low as 
10°K the only rotations possible in ice are about 
the symmetry axis of the molecule. Under these 
conditions the wave functions and eigenvalues 
are no longer those for a spacial rotator with 
three moments of inertia, but correspond more 
nearly to those for a plane (two-dimensional) 
rotator.” The lowest antisymmetric rotational 
state has J=1 and is doubly, not triply, degen- 
erate in the rotation, corresponding to their 
value 3R/4-In 2=1.03 cal./deg. for the zero- 
point entropy which is to be added to the third 
law value. 


19 Giauque and Ashley, Phys. Rev. 43, 81 (1933). 
20 Pauling, Phys. Rev. 36, 430 (1930). 
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1. Heats of mixing for mixtures of primary m-alcohols and 
n-hydrocarbons have been determined down to practically 
infinite alcohol dilution over a temperature range from 
10°C to 45°C. It has been found that (1) the heat of 
mixing is always negative and that the heat absorption, 
per mol of alcohol, increases with increasing dilution of the 
alcohols to a limit of 5800 cal. per mol; that the curves of 
the molar heat absorption against molar alcohol concen- 
tration are identical for all investigated systems—the first 
four primary m-alcohols in n-hexane and n-heptane. 
Furthermore, reports in the literature that n-hydrocarbons 
may be mixed with each other without any heat effect were 
confirmed. 

2. The present experiments and evidence gathered from 
the literature indicate that (1) primary n-alcohols, in the 
pure liquid state, are completely associated to double 
molecules and that when mixed with a hydrocarbon they 
dissociate into single molecules; (2) that the association 
takes place between the two hydroxyl groups of two single 
molecules on account of their electrostatic dipole forces; (3) 
that the observed heat absorption is entirely due to the 
dissociation of double molecules, i.e., breaking up of the 
(electrostatic) hydroxyl bond and that there is no thermal 
effect of the aliphatic groups in the alcohol molecules upon 
the hydrocarbon solvent; and (4) that the equilibria 


between single and double molecules are identical at equal 
molar concentrations for all primary m-alcohols in all 
n-hydrocarbons. 

3. This theory leads to an equation for the partial vapor 
pressures of n-alcohols in n-hydrocarbons: 


log p=log f—(A—Q)/RT+B—, 


where p=the partial vapor pressure at the mol fraction f, 
=the heat of vaporization of pure alcohol, Q=heat of 
mixing at mdl fraction f and temperature 7, and b=a 
universal function of Q. B is determined by the ordinary 
(simplified) vapor pressure equation 


log P= —\/RT+B 


and is constant over the whole temperature range in- 
vestigated. \ and B depend upon the nature of the alcohol, 
while Q and 6 depend upon molar concentration and 
temperature. 

4. It is shown that after b has been determined once as a 
function of Q in a specific case, one can calculate partial 
vapor pressure curves of other systems. 

5. An equation has been developed which permits the 
calculation of the, degree of dissociation and therefore Q, 
for any temperature and concentration. 





HILE extensive experimental material is 
available on the physical properties of 

liquid mixtures, it usually appears impossible to 
obtain an accurate idea of the changes which the 
components undergo during the process of mix- 
ing. Attempts to introduce chemical conceptions, 
such as the formation of double, triple, etc., 
molecules, have frequently been made in the past 
but usually have met with failure; on the other 
hand, no other general method of treatment has 
been offered thus far. However, consideration of 
the experimental data presented in this paper on 
mixtures of primary normal alcohols with normal 
paraffin hydrocarbons appears to justify the 
adoption of a specific and simple picture in this 
particular case. It may be assumed that similar 
simple conditions may be found in other systems 


and further, that the simplicity of the present 
case may allow an approach to the fundamental 
laws governing equilibria in liquid systems. 
Primary normal alcohols, in the liquid state, 
show x-ray diagrams which indicate the existence 
of cybotactic groups containing two sets of 
planes:! one set of planes is produced by the ar- 
rangement of the molecules as double molecules, 
two OH groups being together and the two mole- 
cules forming a straight chain, the other by the 
arrangement of two such double molecules lying 
side by side with the two pairs of OH groups not 
directly adjacent but slipped by one another so 
that the set of planes containing them is not 


normal to the length of the chain. In Fig. 1, the 


‘Stewart, Chem. Rev. 6, 483 (1929). Compare also 
further references there. 





GUENTHER von ELBE 











Fic. 1. Cybotactic group of n-propyl alcohol. 


cybotactic group of n-propyl alcohol is shown, 
from Stewart’s data on the distance of the planes. 
The value of 12.8A for the length of the double 
molecule is estimated from his formula for the 
length of n-paraffin molecules and additional in- 
formation on the dimensions of the OH group.’ 
The distance between the two adjacent OH 
groups is assumed arbitrarily to be that of the 
closest possible approach. 

The head-to-head association is certainly 
caused by the strong electrostatic forces arising 
from the polar OH groups. The dipole moment 
of the OH group itself amounts to about 1.6 
X10-!8 e.s.u.,? which is rather large compared to 
other groups. On the other hand, the separation 
of the poles must be rather small, since the dis- 
tance of the nuclei of the atoms amounts to only 
0.97A.2 It follows that the poles themselves are 
extraordinarily strong and that two molecules 
are linked together with considerable force. In- 
formation on the strength of this electrostatic 
bond may be obtained from a paper by M. 
Dunkel,* who shows that the increase of the heat 
of evaporation of an alcohol over that of the 
corresponding paraffin is fairly constant through- 
out the homologous series and amounts, on the 
average, to 7250 cal. per mol. Since there is no 
association of alcohol molecules in the gas phase, 
this value may be interpreted to represent the 
heat of formation of double molecules. The large 
value of 7250 cal. per mol—comparable to a 
chemical bond—suggests that association in the 
liquid state is practically complete. 


2 Hengstenberg and Mark, Naturwiss. 20, 539 (1932). 
3’ Eucken and Meyer, Phys. Zeits. 30, 397 (1929). 
4M. Dunkel, Zeits. f. physik. Chemie 138, 42 (1928). 
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Coulomb’s law for the potential energy of two 
dipoles, placed parallel with opposite poles ad- 
jacent, is 


V = —2e?(1/x—1/(x?+")!), 


where ¢ is the pole strength, y the diameter of the 
polar group and «x the interpolar distance. If one 
identifies y with 0.97A, which is probably too 
large a value, one may calculate e from the above 
dipole moment and, correspondingly, the poten- 
tial energy as a function of x, as shown in Fig. 2. 
At very close approach, repelling forces will ap- 
pear. If one places the minimum of potential 
energy at X=1A and — V=7250 cal. per mol, 
the dotted curve may represent approximately 
the combined influence of attraction and repul- 
sion. 

Normal hydrocarbons tend to arrange them- 
selves parallel to each other at a distance of 4.6A. 
There is practically no heat effect if one mixes 
two different m-hydrocarbons with each other, 
and the partial vapor pressures (p) above the 
mixture follow Raoult’s law: 


p=fP, (1) 


where f =the mol fraction and P the vapor pres- 
sure of the pure substance. As shown by thermo- 
dynamics, Raoult’s law can hold only in the ab- 
sence of any heat of mixing. 

Dunkel‘ has shown that chemical groups, like 
—CHe, —CHs, preserve their contributions to the 
intra-molecular force in equal strength in differ- 
ent compounds. The aliphatic radical of an alco- 
hol molecule, therefore, has the same effect upon a 
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neighboring molecule as if it were a part of a 
hydrocarbon molecule instead. Therefore, it con- 
tributes nothing to the heat of mixing of n-alco- 
hols and n-hydrocarbons, which appears to be 
entirely due to the dissociation of the hydroxyl 
groups. 

In view of the above considerations certain 
predictions can be made about the effects to be 
observed in the mixing process. At very low alco- 
hol concentrations, the hydroxyl dissociation 
should be complete. Hence, if association of all 
alcohols in the pure state is complete, the heat 
absorbed during mixing should approach a final 
value identical for all m-alcohols in all m-hydro- 
carbons. The present experiments show this to be 
the case and the value found is 5800 cal. per mol, 
which is lower than the value of 7250 cal. quoted 
above because of the higher dielectric constant 
of the liquid medium. At higher concentrations, 
equilibria between double and single molecules 
are established. The degree of dissociation deter- 
mines the molar heat effect. For equal molar 
concentrations, this should also be identical for 
all n-alcohols in all n-hydrocarbons, since the 
concentrations of OH groups are identical and 
the hydrocarbon groups do not influence each 
other. ; 


EXPERIMENTAL PROCEDURE AND RESULTS 


The hydrocarbons used were Eastman “‘prac- 
tical’’ hexane and heptane. They were redistilled, 
especially to insure complete elimination of mois- 
ture. The z-alcohols used were anhydrous 
methyl-, ethyl-, 2-propyl- and -butyl alcohol. 

The liquids were mixed in a Dewar flask, con- 
stantly stirred, and temperature readings were 
taken with a Beckmann thermometer. The alco- 
hol sample was contained in’a special pipette 
which was completely immersed in a weighed 
amount of liquid hydrocarbon and which could 
be opened from the outside. Figs. 3a and 3b show 
two different constructions of such pipettes. In 
Fig. 3a, the alcohol sample is enclosed by mer- 
cury and may be pushed out by turning the 
threaded shaft with the piston on the end. In the 
construction shown in Fig. 3b, two valves, one on 
the top, the other on the bottom, may be opened 
simultaneously, thus releasing the enclosed 
liquid. The amount of alcohol was determined by 
weighing the pipettes before and after filling. The 


















































(a) (b) 


Fic. 3. Sample containers (pipettes) for measuring heats 
of mixing. 


measurements were made by determining the 
temperature drift before and after mixing. Care 
was taken to eliminate heat effects due to evapo- 
ration by keeping the gas space above the liquid 
sufficiently small. The heat capacity of vessel and 
liquid was determined by passing a known cur- 
rent through an immersed resistance wire and 
determining the temperature change and heating 
time. 

In Fig. 4, the molar heats of mixing the above- 
named alcohols in heptane and hexane at 29°C 
are plotted against the concentration (mol per 
cc). The results of experiments with all the 
materials fall on one curve, as expected, and ap- 
proach at infinite dilution a value of Q,,=5800 
cal. This holds for other temperatures as well, the 
only difference being a change in the slope of the 
curves, corresponding to more dissociation at 
higher temperature and vice versa, for any given 
concentration. 


THE DIssOcIATION EQUILIBRIUM IN THE LIQUID 
PHASE 


The equilibrium between double and single 
molecules, at any temperature JT and mol frac- 





GUENTHER 


von ELBE 














~ IN HEPTANE 








°-CH;OH 
*=CzH,OH 


« @ ~IN HEXANE 





# «-C3H,0OH 
9- C4HgOH 





CALORIES PER MOL OF ALCOHOL 


























° 106 y 1o°4 


10> 


MOLS OF ALCOHOL PER CC OF SOLUTION 


Fic. 4. Heats of mixing of primary normal alcohols in hexane and heptane. 


tion f, can be measured by the heat of mixing. 
The degree of dissociation may be given by 


a=0/0. Q,=5800 cal. (2) 


if the assumption is made that no heat is required 
to separate double molecules one from another. 


Since a double molecule resembles a normal 
hydrocarbon and these show no evidence of a 
heat of mixing, this assumption is probably 
justified. 

The results of investigations on the dielectric 
properties of mixtures of alcohols and hydro- 
carbons® may be interpreted from the viewpoint 
developed thus far. The molar polarization P of 
a substance in which the molecules are free to 
assume a perfectly random orientation is repre- 
sented by the Debye equation, written in the 
following form: 


Furthermore, 
[(m?—1)/(m?+2)](M/d) =(41/3)N, (A) 


where =the index of refraction for light of 
infinite wave-length, E=the dielectric constant, 
M=the molecular weight, d=the density, NV 
=Avogadro’s number and yu=the dipole mo- 


5 Smyth and Stoops, J. Am. Chem. Soc. 51, 3312 (1929). 
ps 


ment. The index of refraction should remain 
about the same for dissociated and undissociated 
alcohol, as the electronic structure does not 
change (no “sharing” of electrons). One may, 
therefore, calculate the factor (417/9k) Ny? in (3) 
from Smyth and Stoops’ data. Values thus ob- 
tained for ethyl alcohol in hexane at high dilu- 
tions and ethyl] alcohol in heptane at higher con- 
centrations are plotted as isotherms against the 
mol fraction fin Fig. 5b. In Fig. 5a, corresponding 
heats of mixing curves are shown. Both types of 
curves, from zero concentration upward, drop 
rapidly in the first part. The reason is the same 
in both cases: the degree of dissociation rapidly 
diminishes. A single alcohol molecule possesses a 
strong dipole moment, while double molecules 
are only slightly polar, because of their sym- 
metry. The subsequent slope of the curves in 
Fig. 5b has been discussed by Malsch.® It may be 
pointed out that alcohol double molecules may 
be deformed by transverse vibrations which 
bulge the central OH groups out and cause slight 
polarity. The effect upon an electric field will be 
increased at higher concentrations when several 
double molecules are coupled together in the same 
cybotactic group. 

In order to investigate the quantitative con- 
nection between degree of dissociation, concen- 


6 Malsch, Phys. Zeits. 33, 383 (1932). 
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tration and temperature, the application of the different temperatures and concentrations, where 

law of mass action may be attempted: 2Q.,. = 11,600 cal. represents the heat of formation 

P _. of one mol of double molecules from 2 mols of 

(a?/(1 —a))¢=const. = K, (5) single molecules. This fact permits the calculation 

where c=number of mols per cc. As is shown by Of values of the heat of mixing for any tempera- 

Table I, the mass law itself does not hold for any ture and concentration from one experimental 
isotherm. Surprisingly, however, van’t Hoff’s isotherm of the heat of mixing. 

equation Since the mass law as expressed by Eq. (5) does 

not hold, one may adjust it empirically by adding 

In Ki1/K2=(2Q./R)((T2—T1)/T2T1) (6) an “activity coefficient.” Introducing the factor 


(1.492)4-#)/a/(1.084)?/G-@ in (5), one obtains 


a? (1.492)(1-@)/a 
TABLE I. A pplication of van't Hoff’s equation to the tempera- ox = const. = K*, (7) 


ture dependency of the dissociation equilibrium. l-—a (1.084) 2¢/(l-@) 


holds between the K values for equal a’s at 




















and 

Tec 30 44.4 30444 3042.4 30 42.4 
mol per cc X 104 68 1.7 1.0 2.5 2.0 5.1 5.0 10.7 
a 863 863.755.755.517 517.293 .293 In K,*/K.*=In K,/Ke 
a%</(1—a) X10°=K 36.8 91 22.857.5 11.0 23 6.0 12.9 

Ki\ T1T2 = (2Q../R)(T2— T;)/ TT». (8) 
(i =) R 11700 12200 11300 11800 

Ko} Ti—T>2 


Mean value= 11750 cal. 
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Tec 9.1 30 8.3 42.4 9.2 42.3 | rs 4 
mol per cc X 104 -212 1.01 578 5.6 2.0 13.5 es | } } | 
a .733 .733 495 .495 .245 245 Z 25 ] i 
ec/(1—a) X10=K 4.24 20.3 2.81 27.1 1.59 10.75 a | \ | 
Ki\ TiT2 iL 20 —_—_\—__—. 
In — R 12300 11900 10400 ~ | _ s. | | 
Kz: J 7T:1—T2 8 | | i 
Mean value= 11500 cal. 5 . 
< | Nh | 
Ge | — | 
TABLE II. Comparison between mass action law Eq. (5) and <4 | ha | 
new Eq. (7). Y 05 | ) | 
| | | | | | 
os | | | 
—_— o- A A iw °o ar a2 0s “eT a5 06 a7 a8 a9 10 
15.8 13.5 14.8 


K* x 105 14.1 
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In Table II, the K values of Table I, for 30°C, 
are compared with K* values. The comparison 
extends over a wide range of dissociation and 
shows that Eq. (7) holds satisfactorily. The con- 
stancy of K* obtained here holds for other iso- 
therms as well, of course, because of Eq. (8). In 
Fig. 6 the activity coefficient is plotted against a. 
For a=1 it becomes zero, for a=0 it becomes 
infinitely large. 


THEORETICAL DISCUSSION 


A general theory on equilibria in liquid phases 
does not exist, except for ionic solutions. The use 
of the simple mass action law has frequently met 
with failure in the past and it is not applicable in 
the present case, even if more than one equi- 
librium is assumed. In order to account for any of 
the measured isotherms by assuming triple-, 
quadruple-, etc., molecules it would be necessary 
to consider a very large number of such additional 
equilibria and it would be impossible to account 
for the simple relationship (8) between the iso- 
therms themselves. 

Furthermore, it appears doubtful whether a 
successful theory may be built upon the assump- 
tion of electrostatic forces between the solute 
molecules, as in the case of ionic solutions. In the 
first place, deviations from the mass action law 
are observed even at high dilution, while, on the 
other hand, the electrostatic dipole potential 
rapidly diminishes with increasing distance. 
Second, one will probably find it difficult to fulfill 
the relationship (8). 

Eq. (7) may be formally derived if one is al- 
lowed to assume that the osmotic pressures of 
double and single molecules in separate solutions 
are proportional to their concentrations and that, 
if both kinds of molecules are in the same solu- 
tion together, their respective osmotic pressures 
are increased for an amount proportional to the 
concentration of the other component. Thus, if 
a; and a2 are proportionality factors and ¢; and ¢ 
the concentrations of single and double molecules, 
respectively, the increase in the osmotic pressure 
of the single molecules due to the presence of 
double molecules amounts to ai¢, and the corre- 
sponding increase for the double molecules, as 
compared with their osmotic pressure in a solu- 
tion which is free of single molecules, amounts to 


Q2C}. 


With these assumptions, one may now calcu- 
late the free energy of transforming 2 mols of 
single molecules into 1 mol of double molecules. 
Consider three solutions: one containing single 
molecules at the standard concentration C,,, 
another containing double molecules at the 
standard concentration C:,, and the third con- 
taining single and double molecules together at 
equilibrium concentrations ¢c; and cz. By means 
of semi-permeable pistons we first compress 2 
mols of single molecules from C;, to ¢1; second, 
transfer them into the equilibrium solution; third, 
remove 1 mol of double molecules from the equi- 
librium solution and fourth, expand the latter 
from equilibrium concentration c, to standard 
concentration C2,. In steps one and four the free 
energy change equals simply 2 RT In C,,/c; and 
RT In ¢2/C2;, respectively. In steps two and three, 
work is done because of the increased osmotic 
pressure of each participant. The work which is 
done in step two equals 2 a,v;¢2, where 2 is the 
volume which contains one mol of single mole- 
cules. Correspondingly in step three, the work 
done equals —d2v2¢1, where v2 is the volume con- 
taining one mol of double molecules. The free 
energy of transforming two mols of single into one 
mol of double molecules is, therefore, given by 


AF=2RT In Cis/€1— 2a v1C2 
— RT In C2,/Co+ devel. 
Introducing v:=1/c,; and ve=1/c and trans- 
posing, 
AF=RT In C1.°/C2,—-RT In cy? /cy 
—2(€2/¢1)a1+(€1/C2)a2. (9) 
Since AF, at the temperature 7;, must be inde- 
pendent of any specific equilibrium concentra- 
tions ¢; and C2, one obtains, for constant tem- 
perature 
In €:?/C2+ (2¢2/¢1)(a1/RT) — (€1/¢2) (a2/RT) 
=const. (10) 


Introducing the degree of dissociation, a, by 
means of 


a=C,/C (11) 


and 
1—a=2c/c, (12) 
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where c is the number of mols C,H2,_,OH per 
unit volume, one obtains 


In (2a?/(1—a))c+((1—a@)/a)(ai/RT) 


—(2a/(1—«a))(a2/RT)=const. (13) 


Since the energy of a fully excited degree of free- 
dom increases in proportion to the temperature 
and a, and az possess the dimensions of an energy, 
it is plausible to assume that for different tem- 
peratures the quotient a,/RT and a2/RT remains 


constant: 
a;/RT =const. =In A},- a2/RT 


=const. =In Ao. 


and 
(14) 
Hence 





at A\Gra@la 
x =const. = K*. 
A 92a! (l—a) 


(15) 


1—a 


Introduction of the numerical value 1.492 for A, 

and 1.084 for Az leads to Eq. (7) which is in 

numerical agreement with the experimental data. 
From Eq. (14) one obtains 


a,/T =0.7 cal./deg. per mol and a2/T 
= 0.184 cal./deg. per mol. 


These values may possess a physical meaning in 
connection with the above assumptions. The 
writer believes, however, that further discussion 
of these assumptions would be premature. 


PARTIAL VAPOR PRESSURES 


A knowledge of the heat of mixing makes it 
possible to deduce an equation for the partial 
vapor pressures of n-alcohols in the n-hydro- 
carbons. For this purpose the validity of Raoult’s 
law may be assumed in a modified form: 


p=fxP’, 


where P’ is a hypothetical vapor pressure which 
would be established if pure alcohol were in the 
same state of dissociation as it is in the solution, 
at the concentration (mol fraction) f. The energy 
difference between dissociated and undissociated 
alcohol may be denoted by Li; and one may 
write?” 


(16) 


Q= —Li—(nnyar./Mare.) Le. (17) 


7 Compare Lewis-Randall, Thermodynamics, p. 90 (1923). 
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The second term represents the change in the 
heat content of the hydrocarbon solvent. The 
hydroxyl group of every dissociated alcohol mole- 
cule induces a momentum in neighboring hydro- 
carbon molecules, with the effect that the heat 
content of the solvent increases in proportion to 
the alcohol dissociation, while that of the alcohol 
solute decreases. Substituting, therefore, 


—L,=aq1 and = —(nhyar./Matc.)L2= —ags, 


one obtains 


—L,=agi1=Q0+a¢2. (18) 


Since the specific heats of dissociated and undis- 
sociated alcohol are different because of the crea- 
tion of new degrees of freedom in the dissociated 
state, there is an entropy difference AS. The free 
energy of changing the alcohol modification char- 
acterized by the vapor pressure P’ to the normal 
state characterized by the vapor pressure P is 
equal to 


/ 


F=RT In pahi— Tas=Q- (TAS—ag2). (19) 


It seems again plausible, from the physical pic- 
ture developed above, that AS is determined by 
the degree of dissociation, i.e., Q, alone, and that 
both AS and g2 are independent of the nature of 
either the primary n-alcohol or the paraffin 
hydrocarbon. Combining (19) and (16) and 
transposing 


In p=In f+ln P+Q/RT 


—(1/R)(AS—age/T). (20) 


For the present limited temperature range, AS 
and qg:/T are fairly constant with temperature. 
Introducing, therefore, 


(1/R)(AS—age/T) = 6, (21) 


the term } appears to be essentially a function 
of Q. 

The vapor pressure P is given by the simplified 
vapor pressure equation, which is very satis- 
factory over a limited temperature range 
(ca. 100°), 


In P= —\/RT+B, (22) 


where the constant B may be determined from 














80 


the boiling point, and data for \, in the present 
case, may be found in the literature. 

Combining (20), (21) and (22) one obtains 
finally 


In p=In f—(A—Q)/RT+B-2), (23) 


where \ and B depend upon the nature of the 
alcohol while Q and b depend essentially upon the 
degree of dissociation, i.e., concentration and 
temperature. 

Measurements are available of the partial 
vapor pressure curves of methyl alcohol in hexane 
at 45°C and of ethyl alcohol in heptane at 30°C.® 
Knowing the values of Q and the corresponding 
values of \ and B for both systems, one may de- 
termine } as a function of Q independently for 
both systems. Fig. 7, which is plotted from the 
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Fic. 7. Identity of Q, 6 curves (Eq. (20)) for two different 
alcohol-hydrocarbon systems. 


numerical data given in Table III, shows that 
these functions are indeed identical, as antici- 
pated. 

This means, of course, that having once deter- 
mined b=f(Q) for one system, Eq. (23) will allow 
one to calculate the partial vapor pressure also 
for other systems of this type. This generalization 
seems to be justified despite the fact that the 
comparison is limited to only two cases because 
of scarcity of experimental data. It may be 
pointed out that the two systems under consider- 
ation consist of different alcohols in different 
hydrocarbons at different temperatures and, 
therefore, provide a very severe test. 

Fig. 7 shows that AS is not linear in Q, i.e., a. 
This seems to indicate that the specific heat of 
the solute alcohol does not increase in proportion 
to the dissociation. No suggestion can be made, 


® Ferguson, J. Phys. Chem. 36, 1123 (1932); Ferguson» 
Freed and Morris, J. Phys. Chem. 37, 87 (1933). 


GUENTHER von ELBE 


TABLE III. Application of Eq. (23) to mixtures of ethyl 
alcokol-heptane at 30°C and methyl alcohol- 
hexane at 45°C. 








\ B for logio in mm Hg 


Methyl! Alcohol 9180 8.83 
Ethyl Alcohol 10000 9.16 


Ethyl Alcohol-Heptane, t= 30°C 
0.009 0.020 0.038 0.059 0.100 0.200 





f 
p(mmHg) 17.5* 40 50 55 59.5 63.7 
Q 


5030 3750 2640 1970 1300 770 
b 2.28 1.35 0.73 0.39 0.11 0.00 


Methyl Alcohol-Hexane, t=45°C 
0.030 0.045 0.053 0.070 0.114 0.300 


f 
p(mmHg) 150 229 260 282 300 ~~ 312 
Q 


4750 4080 3780 3280 2270 729 
b 2.09 1.62 1.43 1.17 0.66 0.00 








* Approximate. 






at present, to account for this except, perhaps 
that it may be connected with the peculiar activ- 
ity coefficient in Eq. (7). 

Eq. (23), however, not only allows one to cal- 
culate partial vapor pressure curves from calori- 
metrically measured heats of mixing but to 
predict, from the slope of these curves, whether 
the components are miscible in all proportions or 
separate into two liquid phases. This is demon- 
strated in Fig. 8. If supersaturation would pre- 
vent the separation into two phases, the partial 
vapor pressure, in the range of supersaturation, 
would go through a maximum and a minimum. 
Methyl alcohol and hexane, below 45°C, are no 
longer miscible in all proportions. The partial 
vapor pressure curve of methy] alcohol calculated 


45°C (FERGUSON, JPHYS.CH 26,1125 ; 
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Fic. 8. Formation of two liquid phases in the system 
methyl alcohol-hexane, indicated by calculated vapor 
pressure curve. 
















INTERNAL EQUILIBRIA AND 


from (23) for 40°C indeed goes through a maxi- 
mum and minimum, as shown by the dotted part. 
The actual extension of the range of immiscibility 
differs somewhat from the range over which the 
dotted part in Fig. 8 extends. These limits are 
very sensitive to the Q values used and the lack 
of coincidence of the theoretical and experimental 
range may be accounted for by changes in Q 
within the limit of error in the heats of mixing 
data. 
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The partial vapor pressure curve of the hydro- 
carbon may, in principle, be determined from the 
alcohol curve by means of the Duhem-Margules 
equation 


= (1 —f)d In Pnyar./df =fd In Patconoi/df, 
which, by combining with Eq. (23), becomes 


1 dQ db 
144(——-—). (25) 
RT df df 


(24) 


d| yar. 
aitiademeeits 
df 
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A critical survey has been made of existing data on the heats, free energies and entropies 
of reactions involving aqueous ions. These data lead to remarkably consistent values for 
the ionic entropies. The entropies of 20 positive ions and 12 negative ions have been tabulated. 





ATIMER and Buffington! in 1926 published 
a list of ionic entropies which included many 
of the positive ions and a few of the more common 
negative ions. These values were based on rather 
meager data and consequently an accuracy of 
only two units could be claimed even in the cases 
where the data were most complete. Since the 
appearance of this initial article, an intensive 
program has been in progress in this laboratory 
with the primary object of establishing these 
values more definitely. It is felt that the work has 
progressed to such an extent as to warrant a 
review of the ionic entropy values now available. 

In the initial paper many of the values were 
based on entropies for the solid salts which had 
been calculated by the approximate equation 
given by Latimer? or obtained by comparison 
with similar salts. This uncertainty has since been 
largely eliminated as a result of the specific heat 
measurements made in this laboratory.’ Errors in 
the heats of solution and activities of salt 
solutions are still large in certain cases, but recent 
work by many investigators has reduced these 
uncertainties also. 

The standard state, as in all previous work, of 
the ionic entropies is the hypothetical one molal 
solution, that is, a one molal solution of the ion 
obeying the perfect solution laws and the ions 
possessing the same partial molal heat content 
that they have at infinite dilution. S°29s.1 of Ht is 


1 Latimer and Buffington, J. Am. Chem. Soc. 48, 2297 


(1926). 
2 Latimer, J. Am. Chem. Soc. 43, 818 (1921). 


3 See references to entropies listed in Table II. 
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taken to be zero. The ionic entropies are of course 
partial molal quantities. 

The calculations and the data used have been 
summarized in Tables I and II. 

We have given in Table III the values for the 
entropies of the aqueous ions which appear to be 
the most probable. In general we have not taken 
the simple average of the various determinations 
summarized in Table I but have given weight to 
the data which appear to be the most reliable. In 
addition, we have included the value for fluoride 
given by Latimer and Buffington,! although the 
complete lack of data on the activity coefficients 
of fluorides renders this value only a rough 
approximation. 

In earlier papers the theoretical significance 
and practical importance of the ionic entropies 
have been pointed out: we wish to call attention 
in this paper merely to the consistency which 
now exists in the values determined by wide 
varieties of methods. It should be noted that the 
checks obtained in Table I are, for almost every 
ion, more numerous than the cases listed under 
that ion. Thus, for example, if we use our 
accepted values for iodide, iodate, ferrous, ferric 
and sulfate, we can then calculate by five 
additional methods the entropy of silver ion. 
The maximum deviation from the mean (which 
includes also the error in the negative ion) is 0.8 
e.u. It is also obvious from the data here pre- 
sented that the third law of thermodynamics 
may be applied to salts of non-magnetic ions and 
their aqueous solutions without serious danger of 
errors due to zero point entropy or transitions at 

very low temperatures. 






















of ent fons: e J } ] ] 
of entropy calculations: AH® and AF° in calories ber mole. AS° in calories per degree per mole, (The entro 


TABLE I. Summary 


py data for the elements and compounds used in these calcula- 


tions are summarized in Table TT ) 
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TABLE II. Summary of entropy values used in the calculation 
of Table I. Values in calories per degree per mole. 


SCHUTZ AND J. F. 


G. HICKS, JR. 


TABLE III. Entropies of aqueous ions at 298.1°K. Values in 
calories per degree per mole. 








Substance S°298.1 Substance S°298.1 





Ba(BrOs)2*H20 68.7% F 49.051 
Ba(NOs)2 51.118 K 16.511 
22.511 
35.6% 
19.952 


34.2% 

36. 1% 

F 36.2% 

CaSO: 2H:O0 442 Rb 17.4 
Cc i 10.05 
ital 12H:0 i 25.654 
le " 


Li 
Mg(OH): 
Hg 
Hg2Cle 








Kasplan, J. Am. Chem. Soc. 50, 1863 (1928) gives a value of 
21,887 for this free energy. 

15 International Critical Tables, corrected to 25°C. 

16 International Critical Tables, using the sum of the 
heats of formation instead of the direct determination. 

17(a) Thomsen, Thermochemistry, p. 52, Longmans 
Green and Co., N. Y., 1908. (b) Butler and Hiscocks, J. 
Chem. Soc. 2554 (1926). We have used a mean of these 
values after correcting for dilution and temperature effects. 

18 Latimer and Ahlberg, J. Am, Chem. Soc. 54, 1903 
(1932). 

19 Latimer and Ahlberg, Zeits. f. physik. Chemie A148, 
468 (1930). 

20 Rossini, Bur. Standards J. Research 6, 847 (1931). 

21 Latimer, Schutz and Hicks, J. Am. Chem. Soc. 55, 971 
(1933). 

2 Kline, J. Am. Chem. Soc. 51, 2093 (1929), gives 5 10-” 
as the solubility product of Mg(OH):. 

*3 Lange and Monheim, Zeits. f. physik. Chemie A150, 
349 (1930). 

24 Pearce and Nelson, J. Am. Chem. Soc. 54, 3544 (1932) 
give y+ =0.576 at saturation, 4.81M KCl. 

25 Landolt Bornstein Tabellen, corrected to 25°C. 

26 (a) Seidell, Solubility of Inorganic and Organic Com- 
pounds, D. Van Nostrand Co., N. Y. (1919). (b) Landolt 
Bornstein, Tabellen, II Erganzungband gives y+ as 0.627 
for 4M KCl. Extrapolating to saturated solution we find 
++0.65. 

27 Wiist and Lange, Zeits. f. physik. Chemie A116, 190 
(1925). 

28Lewis and Randall, Thermodynamics and the Free 
Energy of Chemical Substances, McGraw-Hill, N. Y. (1923). 

29 Fromherz, Zeits. f. physik. Chemie A153, 376 (1931). 

30 Landolt Bornstein Tabellen. 

31 Noyes and Beam, J. Am. Chem. Soc. 34, 1016 (1912). 

3% Latimer and Greensfelder, J. Am. Chem. Soc. 50, 2202 
(1928). 
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Ti+ 28.6+0.3 Cl- 13.5+0.1 
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Mgt — 35.5 NO;- 36.9 
ClO, 43.1 
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BrO;— 37.7+41.0 
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33 These data have been obtained from various sources 
including the International Critical Tables and Landolt 
Bornstein Tabellen. Special reference should also be made 
to the work of Richards and co-workers. 

34 Latimer, Schutz and Hicks, J. Am. Chem. Soc. (in 
press). 

3 Greensfelder and Latimer, J. Am. Chem. Soc. 50, 3286 
(1928). These values have been recalculated and a large 
error in the graphical integration of the entropy of 
Ba(BrOs3)2-H:2O has been corrected. The correct value is 
68.7. 

3% Ahlberg and Latimer, J. Am. Chem. Soc. (in press). 

37 Greensfelder and Latimer, J. Am. Chem. Soc. 53, 3813 
(1931). The AH value given in this paper has been re- 
calculated from the solubility data. 

38 Latimer and Ahlberg, J. Am. Chem. Soc. 52, 549 
(1930). 

39 Latimer, Schutz and Hicks, J. Am. Chem. Soc. 55, 971 
(1933). 

40 Backstrém, J. Am. Chem. Soc. 47, 2432 (1925). 

4t Latimer, Hicks and Schutz, J. Chem. Phys. 1, 424 
(1933). 

Latimer, Hicks and Schutz, J. Chem. Phys. 1, 620 
(1933). 

43 Latimer and Hoenshel, J. Am. Chem. Soc. 48, 19 
(1926). 

44 Giauque and Overstreet, J. Am. Chem. Soc. 54, 1741 
(1932). 

4 Giauque, J. Am. Chem. Soc. 52, 4816 (1930). 


46 Giauque and Wiebe, J. Am. Chem. Soc. 50, 101 (1928). ‘ 


47 Giauque and Wiebe, J. Am. Chem. Soc. 51, 1441 
(1929). 

48 Giauque and Ashley, Phys. Rev. 43, 1 (1933). 

49 Kelley, U. S. Bureau of Mines Bull. 350 (1932). 

50 Giauque and Archibald, unpublished data. 

51 Giauque and Johnston, J. Am. Chem. Soc. 51, 2300 
(1929). 

82 Nernst, Ann. d. Physik 36, 395 (1911). 

53 Giauque and Meads, unpublished data. 

54 Eastman and Milner, J. Chem. Phys. 1, 444 (1933). 
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On the Statistical Basis of the Theory of Electrolytes 


Otto HALPERN, Physics Department, New York University, University Heights 
(Received December 14, 1933) 





In the following paper the author investigates how far the _ these limitations has been responsible for certain paradoxes 
various statistical-thermodynamic methods hitherto pro- _ in the literature on this subject. A new fluctuation formula 
posed are suitable for application to strong electrolytes. for the electric potential of the ions is derived and con- 
The significance of the ionic radius is, in many instances, clusions are drawn from it as to the correctness of existing 
far more important than has hitherto been assumed. When _ theories. In connection with this last a deduction proposed 
the influence of the ionic radius is negligible conditions are by Fowler is reviewed. Kramers’ theory of electrolytes is 
derived from dimensional considerations which con- tested with respect to its agreement with the theorems here 
siderably restrict the functional form of the free energy. developed. Finally, a review of the statistical and thermo- 
Furthermore, proceeding from conditions of integrability, | dynamical basis of Bjerrum’s theory of ionic association is 
thermodynamics appreciably limits the number of possible _ presented. 

expressions for the ionic potential. The disregarding of 





1. INTRODUCTORY REMARKS Kramers? which proceeded from the rigorous 
NEXPECTED progress has been made in phase integral of Gibbs’ canonical ensemble and 
the theory of strong electrolytes in the from this aimed at obtaining Debye’s formulas. 
course of the decade since Debye’s theory was The second group contains first and foremost an 
first put forward. Whereas (assuming the cor- important inquiry by Gronwall, La Mer and 
rectness of Debye’s deduction for the mean ionic Sandved* who proceed from the same differential 
potential) the thermodynamic part was prac- equation as Debye but aim at putting it from 
tically disposed of, Debye’s conception of the the start on a stricter mathematical basis by cal- 
ionic cloud, on the other hand, gave rise to the Culating the higher approximations, a procedure 
discovery and interpretation of a number of which would be unfeasible on principle if the 
effects in the field of conductivity. The knowledge ionic radius were disregarded. Finally, the third 
thus gained must be regarded as an addition of 8'OUup is composed of an investigation by Bjer- 
permanent value to this field. rum‘ and of subsequent work following the 
The basis for this development is provided by Program outlined in it. In this were introduced, 
Debye’s method of calculating the mean ionic by way of completing Debye’s mode of approach, 
potential. That this method does not follow di- the conceptions of ionic association which are to 
rectly from the recognized laws of statistics but be explained below. 
must, on the contrary, first be justified by means The result of these investigations as given by 
of statistical principles, has already been empha- the authors themselves, can in its main features 
sized by so many authors that it is hardly neces- be summed up as follows: A concentration can 
sary to go further into this point at present. The @!ways be found below which the Debye theory 
matter lies about as follows: because of the re- ‘etains its validity as a limiting law. For higher 
peatedly observed agreement in the field of dilute Concentrations deviations from Debye’s limiting 
aqueous solutions between Debye’s theory and aw enter which vary according to the different 
experiment a great variety of methods has been authors. In this connection the Debye limiting 
tried to vindicate this agreement in theory. The aw is defined by the statement that the excess 
efforts made in this direction are essentially divis- 2H. A. Kramers, Proc. Amst. Akad. Sci. 30, 145 (1927). 
ible into three groups. To the first of these belong 3 T. H. Gronwall, V. K. La Mer and K. Sandved, Phys. 
the investigations carried out by Fowler! and  Zeits. 29, 358 (1928). 


4N. Bjerrum, Copenhagen Akad. Inst. fys. Medd. 7, 9 
1R. H. Fowler, Statistical Mechanics, Chapters 8 and 13. (1926). 
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free energy of an ionic solution is given by the 


formula 
F=-); Nie2«/3D, (1) 


k= (4r>0:Nie?/DRTV)?. 


The symbols entering into this formula require 
no further explanation. 

It is the purpose of the present investigation to 
test in the first place the correctness of these re- 
sults, secondly to show how far it is possible to 
draw from general properties (conditions of 
dimension and integrability, etc.) conclusions 
which must be valid for all theories without ex- 
ception, and finally to test the validity of certain 
corollaries of the theory with the aid of general 
thermodynamic conditions. 


2. DIMENSIONAL CONSIDERATIONS AND THERMO- 
DYNAMIC CONDITIONS OF INTEGRABILITY 


We shall consider a solution composed of 
various ions of charge ¢«;. The number of ions of 
the 1-th group will be denoted by Nj, the total 
number of all ions by N, the total volume of the 
solution by V, the dielectric constant by D and 
the radii of the ions of the i-th kind by a;. The 
difficulty inherent in the use of the macroscopic 
dielectric constant for the treatment of the 
microscopic forces between the ions has already 
been emphasized by several authors and will not 
be gone into here, since we are merely able to 
recognize its existence but not to overcome it. 

Between the ions act the Coulomb forces and 
repelling forces of the van der Waals type, which 
prevent the uniting of two opposite charges. We 
have already partially taken these forces into 
account by the introduction of ionic radii. The 
electrical potential energy of two ions of the 7-th 
and k-th kinds is then 


€i€k, Drix, 


in which 7;, denotes the instantaneous distance 
between the two ions. For the purpose of the fol- 
lowing considerations we shall assume the ions 
to be consecutively numbered from 1 to N, and 
the 7-th ion to have the charge e;. There are then, 
accordingly, groups of N; ions present which all 
have one and the same charge ¢;, although it is 
not necessarily to be assumed that any two ions 
in the solution have the same charge. The total 
-excess free energy of the ionic solution, insofar as 
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it is due to the forces between the ions, is then 
known to be given by Gibbs’ phase integral 


e~(FIKT) vem f dey »-dzye—F/kT, (2) 


In E both the Coulomb and the van der Waals 
forces are to be included. If the potential energy 
of the van der Waals forces (short range forces) 
be put down as Ey, then E equals 


E=Ect+Ew=Ewt3 ew €:€%/ Dr ix. (3) 


In this expression for the potential energy the 
proper electric energy ¢,;?/2Da; of each ion is 
omitted as being an uninteresting constant. If 
E be written in the form 


E=Ewt} Dewi, (4) 
then 
Vi= Die €x/Drir (5) 


denotes the instantaneous potential produced by 
all other ions at the position of the 7-th ion. 
In any electrostatic system the relation® 


¥i=0E¢/d¢; (6) 


which follows directly from (3) and (5) is valid. 
The mean value of the potential ¥; then results 
on taking the statistical mean value of the ex- 
pression dE ¢/de;. This statistical mean value is 
obtained by differentiating both sides of (2) with 
respect to ¢e;. Thus 

OF JS dx;: ° -dzy (OE c/d€;)e"-2w-2 o)/k? 


i= 
O€; V¥ 





(7) 
From this follows at once by a second differen- 
tiation the fundamental relation 


OY :/de, = Ii /de:. (8)* 


In the absence of van der Waals forces (vanish- 
ing ionic radius) the phase integral clearly di- 
verges, since all ions of opposite sign simply 
neutralize each other. A finite ionic radius (pres- 
ence of van der Waals forces) is thus of decisive 


5 The method employed at this point and subsequently is 
closely connected with that used in an earlier investigation. 
Cf. Gross and Halpern, Phys. Zeits. 26, 403 (1925). 

* All differentiations with respect to ¢; are, of course, 
carried out at constant temperature; the resulting rela- 
tions for y therefore refer to isothermal properties of the 
ionic potentials. 
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importance for the statistical problem. There is 
therefore no concentration at which the value of the 
free energy is independent of the ionic radius (1.e., 
of the magnitude of the van der Waals forces). It 
is merely possible, as shown by experiment, to find 
cases 1n which, for a given range of values of the 
ionic radius, the free energy is not appreciably 
dependent on the actual value of this ionic radius. 
In such cases we can write, as an approximation 


E=Ec=} Liew (3a) 
and thus obtain by means of (7) the relation 
E=F-kT(dF/a(kT)) =3 Di (OF /deiex (9) 
or, by a simple transformation (n;=«,), 
Pati +- fe. (10) 


d(kT) an i 


The solution of this Eulerian partial differential 
equation is 
F=kTG(ni/kT), (11) 


in which G signifies an arbitrary function of the N 
arguments 7:/kT. G must have the dimension 0; 
since 9;/kT is of the dimension of a length, G can 
only be made dimensionless by the introduction 
of the factor V-}, so that Eq. (11) takes the form 


F=kTG(ni/kTV'). (12) 


This dimensional equation previously has been 
obtained by O. Klein® using a statistical method. 

These dimensional considerations may, how- 
ever, be carried considerably further: As is read- 
ily seen from (3), E is a symmetrical function of 
the «;, This symmetry must accordingly also be 
valid for F. We may therefore add to (12) the 
further condition that G depends symmetrically 
upon the arguments 7;. 

Further, as long as we ignore surface phenom- 
ena, the free energy, as always in thermodynam- 
ics, must be proportional to the total number of 
the ions, so that the N; enter in the expression for 


6O. Klein, Medd. fr. K. Vet. Nobelinst. 5 (1919). 
A derivation of this same equation given by P. van 
Rysselberghe, J. Chem. Phys. 1, 205 (1933), with the aid of 
the virial theorem can not be regarded as satisfactory, 
since in the derivation (Cf., e.g., the transition from (12) to 
(14) and from (21) to (23)) quantities are treated as 
independent where as they are actually dependent on the 
variable of differentiation. 
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F only in the form N;f(N;/V). For the following 
considerations it is of decisive importance to ob- 
serve that F does not contain the number JN; of 
ions of the i-th class explicitly but only through 
the fact that precisely N; ions have all the same 
charge «;. We thus obtain a further limitation, 
since F has to be of the form 


UE 
kTV} 





Ni 
F=) ;—H; 
Ee 


H;=H;, 


N; 
=rWak(—), (13) 
V 


(on account of symmetry !). 


This limitation might appear at first sight to be 
of subsidiary importance, but this is by no means 
the case.” 

It must once more be emphasized that all con- 
siderations which disregard the ionic radius can 
not be of general validity but are only justified 
when it is known that the ionic radius lies within 
a certain range and that the results are practically 
independent of its actual value. 


3. APPLICATIONS OF THE RELATIONS (8) AND (13) 


Let it be first assumed, as is done in the in- 
vestigations of a number of authors, that F is 
only a function of the ionic strength. This is de- 
fined by the expression 


F=)): Nie?f(LiNie?), 


>: Nie? =ionic strength. 


(14) 


This assumption is permissible, since according 


7One might cast doubt upon the correctness of this 
argument by referring to the fact that perhaps the number 
N enters into the expression for the free energy not only by 
the number of the charge carriers but also through the 
volume V which is available to all the particles. This, 
however, is not the case: Let us assume that the N; enter in 
some way into the expression for F not only through the e; 
but also through V. Let us then compare the excess free 
energies of two solutions, one containing N ions and the 
other M additional ions carrying infinitesimally small 
charges (in the limiting case the charge 0). The excess free 
energies must then obviously be equal since the ions with 
the charge 0 can have no effect. They would be unequal 
however if V were to contain the number of ions as a 
factor. It must be borne in mind that we are here dealing 
only with excess free energies as defined above. Failure to 
notice this fact is responsible for a conflicting result 
obtained by Kramers (cf. §5 of this note) which therefore 
must be held to be invalid. 
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to it F becomes a symmetrical function of ¢,?. In 
order to give (14) the form (13) we write 


LiN ie; DiNie? 
" a ). 
yi kTV? 


F 





(15) 


It will be seen at once that (15) satisfies the condi- 
tion (13), only when F is of the form® 


F= CY Nie?/D) (>; Nie?/DkRTV) 4, (16) 
C=pure number. 


From general dimensional considerations there- 
fore we obtain the following result: Assuming that 
the free energy is a function merely of the tonic 
strength and assuming further that the tonic radius 
hes in a range in which its influence may be disre- 
garded, then a limiting law for the free energy 1s 
always valid and this law can differ from that of 
Debye only by a numerical factor. 

A similar conclusion can be drawn in a related 
case. Suppose it is known, for example, that the 
mean potential ¥; is of the form ¥;=e:f;, in which 
f; denotes any symmetrical function of the 7; 
(possibly with the omission of the variable 7;). 
The condition (8) then leads to the relation 


€:0f ;/ Oe, = €x0f ./O€:. (17) 


This can only hold true when f; contains the e; in 
the combination 


fi=fe=f(Deer). (18) 
For, in this case, 
2eiexf’ = 2exeif’, (19) 


in which the prime denotes differentiation with 
respect to the argument (> e;?). From this point 
on the argument runs as above. 

A considerable number of proposed theoretical 
formulae do not satisfy the general conditions 
here given or do so only in special cases obtained 
by artificial restrictions on the parameters, 
_restrictions for which no physical reasons what- 
soever can be assigned. Thus, for example, the 
expression previously given by Debye for the 
mean potential of an ion, as a function of the 
ionic radius, 


¥i= —ex/D(1+xa,), 


8 It follows from (2) and (3) that E and therefore F 
contain the ¢ only in the combination ¢/D. 
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is only in agreement with (8) if a;=a; be valid’ 
The same argument applies to the more con- 
sistent method of integration of Debye’s differ- 
ential equation given by Gronwall, La Mer and 
Sandved, although it becomes clear, when the 
higher approximations are taken into considera- 
tion, that even in the case of equal radii the condi- 
tion (8) is in general not fulfilled, wnless all the 
ions carry the same charge. Thus, for example, 
Gronwall, La Mer and Sandved obtain for the 
potential of an ion originating from a uni-divalent 
electrolyte the expression 


@ 


Vi=Lim(ePm1/as") Gms (20) 
in which ¢,,(e7) denotes a symmetrical function 
which is of no importance for our further deduc- 
tions. It will be seen at once that (20) no longer 
satisfies the condition (8) in higher approxima- 
tions, since if we write ¥; in the form 


— @ 
vi= Lav i™ 
1 


the relation 
Oy &™ /de; = OY, ™ /de; 


is no longer valid for m>1. In this connection it 
must be emphasized that the calculation ac- 
cording to Gronwall, La Mer and Sandved is 
more consistent and that it alone gives a justifi- 
cation of Debye’s method of calculation. The 
fact that the thermodynamic integrability condi- 
tions can not be satisfied points to serious in- 
herent weaknesses of the theory which only ac- 
cidentally do not become manifest in the first 
approximation (with ionic radii equal). All 
theoretical attempts to treat the case of unequal 
ionic radii on the basis of the Debye theory, have, 
so far as the author knows, led to formulas which 
fail to satisfy the integrability condition (8). 

For practical purposes this limitation has been 
of no great importance, since ionic radii have been 
introduced, at least in aqueous solutions, only 
for the purpose of correction, and differences 
between individual ionic radii lie outside the 
range of accuracy of the theory. But it must be 
observed, in order to appreciate the argument, 
that the equality of the ionic radii enforced by 
thermodynamics is in fact a physically absurd 
limitation. 
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A paradox emphasized in the paper by Gron- 
wall, La Mer and Sandved® quoted above also 
finds its explanation in this manner. If we wish 
to calculate the free energy of an ionic solution 
from the given potentials yi, it is necessary on ac- 
count of (7) to carry out the integration of the 
following system of N partial differential equa- 
tions for the one unknown quantity F 


OF /de=Vi 
(t=1, 2---N). 


(21) 


Since for e¢;=0 the free energy is equal to 0, we 
can in every case take e;=0 as the lower limit of 
integration. The upper limit of integration is 
given of course by e;=«;. Since we are here deal- 
ing with N independent variables it is possible to 
carry out the integration in any number of ways. 
The Eggs. (21) and (8) invariably lead to the same 
result, since it is clear that (8) expresses the nec- 
essary and sufficient conditions for the inte- 
grability of (21). If (8) is not fulfilled then this is 
no longer the case. Gronwall, La Mer and Sand- 
ved carried out the integration by causing (at 
least for symmetrical electrolytes of uniform 
ionic radius) all ionic charges to increase at an 
equal rate so that during the process of integra- 
tion all ions had the same charge at any one time. 
These authors thus maintained, during the inte- 
gration process, equality of all charges concerned, 
a condition on which the fulfillment of (8) had 
been shown to depend. Their result differs from 
that obtained by H. Mueller, since the latter 
carried out a different method of integration. He 
began by causing the charge of the k-th ion to 
increase from 0 to ex, while the charges of the 
other ions were kept constant. During this proc- 
ess there were thus present in the solution ions 
“of different valencies,” so that the integrability 
conditions (8) were not fulfilled; the deviation of 
his result from that of Gronwall, La Mer and 
Sandved, in whose work the integrability condi- 
tion was (artificially) fulfilled, is thus readily ex- 
plained.'° 


* Gronwall, La Mer and Sandved, Phys. Zeits. 29, 367 
(1928). 

10 The significance of the integrability condition and its 
use to explain the charge paradox has already formed the 
subject of a communication previously made, as a pre- 
liminary statement, by the present author. Cf. Phys. Rev. 
43, 495 (1933). 





4. FLUCTUATIONS OF POTENTIAL 


Debye’s methods of calculation, which have 
proved successful in so many instances, consist, 
as is well known, in determining the distribution 
of the charge in the vicinity of a particular ion by 
combining the Poisson equation with a somewhat 
free analog of the e-theorem of Boltzmann. 


AY = —exp/D=—(en/D)X poet”, 


This use of the Boltzmann theorem urgently re- 
quires justification by strict statistical principles. 
If the potential, which also occurs as argument of 
the exponential function, were determined solely 
by the particular ion, then this application would 
be justified. In the case with which we are here 
concerned it is precisely the contributions of the 
other ions to this potential which are of decisive 
importance and for these the simple formulation 
of the Boltzmann theorem here introduced is not 
valid. These objections were already raised some 
time ago by a number of authors and the investi- 
gations of Kramers and Fowler, mentioned in the 
introduction to this paper, aimed primarily at 
dealing therewith. 

As to the magnitude of the error introduced by 
the use of the Boltzmann theorem, no theoretical 
statement has hitherto been possible although 
some information on this point can be obtained 
with the aid of astrictly valid fluctuation formula 
for the potential. For this purpose we proceed 
from the relations (2), (4), (6) and (7) which are 
strictly valid even in the case of the presence of 
van der Waals forces, and obtain by differentia- 
tion of (7) with respect to e;. the relation 


OY: a (dEc 
v*% : = f aes --dzy —|- e(F-Ew-Ec)/kT | 
€i €; de; 
(23) 


(22) 








The expression on the right, however, becomes 
simplified in view of (Cf. (6) and (7)) 


CEc/ de,” = 0, OF, 0€; = Vi 
to 


dvi 
O€; 
X (Wivi— V2) = V9 (1/RT) (Vi — V2). 


In this connection, as is once more emphasized, 
¥; is merely the potential produced by the other 





Vy 


1 
=— | dey . -dzye(F—-E2W-EcyikT 
kT- 


(24) 
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ions at the position of the i-th ion. For the quanti- 
tative evaluation of this fluctuation formula it is 
necessary to know the expression for d¥;/de;. 
For our purposes it is perfectly sufficient to re- 
gard it as empirically given, by osmotic or activ- 
ity measurements. Then, as a close approxima- 
tion, we have in dilute solutions 


dY:/des=Yi/es 
and therefore 
VeVi = — (kT Vi. 
Now Debye’s theory has hitherto been applied 


to aqueous solutions practically only for such 
ranges in which 


(25) 


ewi/kT&1. (26) 


Indeed the estimates to be discussed later, as to 
the range of validity of the Debye formulas, 
furnished as a necessary limitation precisely the 
relation (26). In this case, however, we see from 
the strictly valid fluctuation formula that the po- 
tential produced by all other ions exhibits 
fluctuations which are very great compared with 
its mean value. Moreover this ratio—and this is 
perhaps the most important point in the present 
analysis—is made not better but worse when we 
pass over to smaller concentrations. For then Vi 
approaches 0 as N/V, but ¥? only as (N/V)}?. 
If possible at all, it is thus only at higher concen- 
trations that the use of the mean value of y in 
(22) can to some extent be justified; but then 
the van der Waals forces must be explicitly taken 
into account. 

The higher fluctuations of ¥;, ¥,’, etc., can be 
determined in an analogous manner; since, how- 
ever, they do not lead to further restricting condi- 
tions, they will not here be explicitly given.* 


5. REMARKS ON THE INVESTIGATIONS 
BY KRAMERS AND FOWLER 


The result arrived at in the preceding para- 
graphs conflicts with a relation deduced by 
Fowler also from statistical mechanics, according 
to which it would appear that at sufficiently low 
concentrations, Debye’s method of combining 


*Similar applications of the fluctuation-formulae can 
be made to the problem of electric and magnetic polariz- 
ability (Clausius-Mosotti relation and Lorentz’s inner 
field). Cf. a preliminary note of the present author, Phys. 
Rev. 40, 134 (1932). 


the Poisson equation and the Boltzmann theorem 
was correct. Fowler! derives the following rela- 
tion: 


emAyitAmpem/D 


-leL) 


The left side of this equation requires no further 
explanation; for the exact meaning of the density 
jp here introduced, reference should be made 
to Fowler’s own arguments. W denotes the 
energy which the m-th ion situated at the point 
x, y, Z assumes, under the influence of the particu- 
lar ion and of all other ions, }>(@W/dx)? the 
square of the gradient of this energy, and W,, 
the mean value of W. In Debye’s theory the 
right side of (27) is made equal to 0. This is 
permissible according to Fowler, provided the 
terms on the right side are small compared with 
those on the left side. He shows by explicit calcu- 
lation that }°zy:(9W»/dx)*1/kT is, in fact, sufh- 
ciently small when 


Eni /RTK1 


} (27) 


holds. This implies that in Fowler’s opinion we 
may, in fact, make the right-hand side equal to 
0, at very low concentrations. 

In this connection, however, it has not been 
taken into consideration that (@W/dx)? in view of 


OW A@Wn\? aW? OWn\? 

et ie Oe Yt 

Ox Ox Ox Ox 
is invariably greater than (dW,,/dx)? and in the 
present case even becomes of a higher order of 
magnitude. This can be demonstrated in the 
following manner. dW/dx denotes the instan- 
taneous force exerted in the x-direction upon a 


charge €» at the point x, y, z. The mean value of 
the expression 


(1/RT €m?) ¥ 2y2(0W/dx)? 


is thus equal to the mean of the square of the 
field strength at the point occupied by the charge. 
The order of magnitude of this quantity (£)? can 
be calculated from the general relation 





A =(DV/8n)(E)?, (28) 


in which A denotes the free electric energy of the 
solution. This free electric energy contains as 
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largest term the proper energy of the ions which 
for spherical ions is approximately given by 


A = > Nie? /2Da. 
We thus obtain 
1 OW\? €m?—— 81€m7A 
—r(—) -2@'-——_ 
kT Ox kT DVkT 
Siem? Nie,” Gn” 


= = K 
DkTV 2Da 2aD 


(29) 


(30) 





The other terms in the free electric energy A, due 
to the interaction of the ions, will be disregarded 
in this connection, since we have already shown in 
the preceding paragraph that they are of an ap- 
preciable order of magnitude. At all events, pro- 
vided the concentrations concerned are not too 
high, the following holds:" 


| Ke mW | = €m2x?/DKA/RT Y ey2(9W/dx)? 
=€m’k?/2aD (31) 


for all xK1/a. This shows again that it is not 
justifiable to make the right side equal to 0 and 
that going to lower concentrations by no means 
improves matters. 

Kramers has attempted to evaluate the Gibbs 
phase integral directly by a very ingenious 
method. Unfortunately the assumptions and 
omissions made in the course of the integration 
can hardly be estimated in their effect on the 
final result. We shall, therefore, not enter into a 
discussion of the various details and shall con- 
tent ourselves with drawing attention to only 
two points of general significance. As the result 
of his calculations, Kramers obtains for the free 
energy an expression of the form 


F=—); Nee?(x/3D+CeP+---). 


As will be seen, the free energy is here only a 
function of « and therefore of the ionic strength; 


(32) 


4 It will be seen from (30) that 2(@ W/dx)? is proportional 
to x? while all the terms on the left side of (27) are pro- 
portional to «*. L. Onsager in a lecture given before the 
American Chemical Society in Washington, 1933, has 
pointed out that the consideration of the van der Waals 
forces would reduce the fluctuation terms. But this could 
not make the right side of (27) of smaller order of magni- 
tude (proportional to x‘, as Fowler supposed). It follows 
from the general fluctuation-formula (25) that at low 
concentrations the fluctuation terms always become 
dominant as compared with the retained average values. 


in spite of this other terms occur in (32) besides 
those giving the limiting law (1). This is contra- 
dictory to the rigorous dimensional considera- 
tions previously given (Cf. in particular footnote 
after (13) and §3), so that for this reason, if for 
no other, the formula for the free energy can not 
be correct. 

The second important point is to be found in 
the way in which the ionic radius (the van der 
Waals forces) is treated. Kramers himself 
emphasizes the fact that the phase integral is 
divergent for infinitely small ionic radius. Never- 
theless he expresses the opinion that it is always 
possible to find a concentration at which the value 
of the free energy is independent of the ionic 
radius. This can not be admitted. Since the inte- 
gral becomes infinite for a =0 there can never be a 
concentration at which it is independent of the 
value of the ionic radius. As explained above, we 
can only hope to make the free energy indepen- 
dent of the particular value selected for the ionic 
radius if this value is contained in a range fixed 
from the outset. Experiment seems to indicate 
that this is possible in the case of aqueous solu- 
tions or other solutions of sufficiently high 
dielectric constant. The universal existence of a 
limiting law can not, however, be admitted. This 
statement may be illustrated by a simple estimate 
obtained from Gronwall, La Mer and Sandved’s 
integration of Debye’s equation. Let us consider 
an electrolyte in a solution of low dielectric con- 
stant (f.i. dioxane D~2). Now Gronwall, La Mer 
and Sandved” give an expression for the activity 
coefficient from which the regions of concentra- 
tion can easily be estimated in which the terms 
containing the ionic radius may be disregarded in 
the expression for the free energy, in comparison 
with the limiting term. The condition for this 
can be reformulated as follows: 


9 


fore) ka 


1 = : 33 
> ——e ee sais 





The sum 





< Ka ( e )" 
1 (2m—1)(2m+2)! \DkTa 


can easily be calculated for sufficiently large 


12 Gronwall, La Mer and Sandved, Phys. Zeits. 29, 376 
(1928), formula 69. 
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values of e/DkTa. An asymptotical expression is 
thus obtained 


1 /DkTa\? 
E=3( ) ee /DkTa_ (34) 


e 





For a reasonable value of the ionic radius 
a~5.108 cm 
there results the condition 
KK (1/2a) (e2/DkTa)%e—2/ Pk 74 (35) 
or for the ionic concentration in mols p. lit. 
cK10-*, 


As will be seen, there can be no question of a 
valid limiting law for any reasonable concentra- 
tions. At the same time, of course, the method 
adopted by Gronwall, La Mer and Sandved 
though strictly consistent in itself is nevertheless 
subject to the objections, given in §3 and §4. For 
such low concentrations only the fluctuation term 
in (25) is decisive, the application of the Boltz- 
mann theorem therefore is even for approximation 
purposes, completely illusory. 

The result of these arguments can perhaps be 
summed up as follows: Contrary to general 
opinion, it is impossible by a reduction of con- 
centration always to reach a region in which the 
limiting law holds. While the Debye theory as 
extended by Gronwall, La Mer and Sandved 
permits to do so, though for media of low dielec- 
tric constant only at concentrations without 
physical interest the consideration of the fluctua- 
tion terms excludes this possibility. There would 
be nothing like a limiting law in water either if the 
ionic dimensions happened to be, let us say, one- 
fiftieth of their actual value. 


6. ByERRUM’s IONIC ASSOCIATION 


In the method proposed by Bjerrum for the 
calculation of the free energy with the aid of 
‘ionic association” the following essential points 
of view stand out. Whereas Debye characterizes 
the interaction between a particular ion and the 
other ions by the introduction of the concept of 
an ionic cloud in which the charge is continuously 
distributed with variable density, Bjerrum at- 
tempts to take account of the atomistic concen- 
tration of the charge, by considering the isolated 


interaction between two ions. He thus subdivides 
the neighborhood of each ion into two parts; in a 
region close to the ion he determines the probabil- 
ity of finding another ion while outside a certain 
radius R Debye’s method is employed. Bjerrum 
himself has already pointed out that this division, 
and the choice of the radius R contains arbitrary 
elements. If a pair of ions is located within the 
radius R, this pair is regarded as a ‘‘molecule”’ 
and is not to be considered as contributing any- 
thing to the ionic atmosphere. 

This method is mathematically completely 
analogous to another which for decades has been 
used to calculate the equation of state of gases in 
the presence of van der Waals forces.'® The 
physical justification of this procedure was given 
by the fact that the intramolecular forces con- 
cerned decreased with a high power of the dis- 
tance. Precisely for this reason the method breaks 
down for Coulomb forces. Bjerrum’s calculation 
of the ionic association also differs in two points, 
to its disadvantage, from the determination of the 
van der Walls correction: The sphere of action of 
two atoms within which they can be regarded as 
united to a “molecule,”’ is actually indeper lent 
of the temperature, whereas according to Bjer- 
rum we find it given by the expression 


R=ee,/2DkT. (36) 


The dependence on T in (36) can not be avoided 
since the dependence on D is quite essential and 
according to statistical mechanics (Cf. (2) and 
(3)) the excess free energy divided by kT is only a 
function of the expression 


€:€,/DkT. 


In Bjerrum’s case this has the troublesome conse- 
quence that the Helmholtz equation, which con- 
nects energy and free energy, is no longer valid. 

To prove this we refer to Bjerrum’s formula for 
the excess free energy f of a pair of ions forming a 
molecule 


R 
e-UIFT) (4.7 R3/3) = f Amr2dre(eierl DET) 
P (37) 


(aKR). 


Differentiation of (37) with respect to T leads to 


13 Cf. R. H. Fowler, Statistical Mechanics, Chapter 9. See 
also further references. 
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+4eR%e-(iet/DETR)(QR/IT) (38) 


or after a slight transformation, denoting by @ 
the average excess energy of the ion pair, to 


f—T(0f/dT) =a4+3kT (1 —el—(eiek/ DR) kT) 


(@R/dT = —R/T). (39) 


The additional term 


3kT(1 — elf—(eiek/DR)]/kT) 


is, in general, not small but of the same order of 
magnitude as 7%. ; 
The second point which leads to difficulties in 
quantitative applications is the question how the 
“proper-volume” of the “‘molecules’”’ is to be 
taken into account. Ions spaced at a distance 
smaller than R are to be regarded as ‘‘associated.”’ 


For solutions of small dielectric constant and 
ol 
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particularly for multivalent ions, R assumes 
very large values, precisely on account of the 
long range of the Coulomb forces. Thus, for 
example, for univalent electrolytes in dioxane, 
at room temperature, R takes the value 1.4 
x 10-* cm. The volume correction in the expres- 
sion for the free energy of the non-associated ions 
which can easily be made in the case of van der 
Waals forces can no longer be carried out in this 
instance.'* It will also be impossible, in general, 
except for infinitesimally small analytic concen- 
trations, to state to which of the other ions a par- 
ticular ion belongs. The difficulty involved in the 
fact that the forces in electrolytes are long range 
forces can not be overcome by the application of a 
method capable of furnishing valid results for 
short range forces only. 

To his colleague, Professor Doermann, the 
author wishes to express his thanks for revision 
of the manuscript. 


144 The argument given here will be applied in a discussion 
of electrical conductivity in media of small dielectric 
constant (Cf. Gross and Halpern, this Journal, to be 
published). 
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Dielectric Constant Studies. III. Aqueous Gelatin Solutions 


A. L. FErGuson, J. GILBERT MALONE AND L. O. CAsE, Chemistry Laboratory, University of Michigan 
(Received December 11, 1933) 


Measurements have been made of the dielectric constant 
of gelatin solutions. A test of the Bjerrum ‘‘Zwitter-ionen” 
theory was furnished by a study of the effect of pH on the 
dielectric constant as measured by the Drude method; the 
data obtained support the theory. Through the application 
of an equation suggested by Marinesco, an attempt was 
made to determine the amount of water fixed by gelatin; if 
this equation is valid for gelatin the results obtained for 


absorbed water must be interpreted as indicating that the 
water bound by the gelatin is still free to show its dielectric 
properties. By means of another equation given by 
Marinesco relating dielectric absorption with molecular 
weight, values were calculated for the molecular weight of 
gelatin; the results obtained are consistent with those given 
by other methods. 





ECAUSE of the probable dipole nature of 
proteins, it should be possible to obtain 
some valuable information concerning their 
properties through dielectric constant measure- 
ments. Thus, a test of Bjerrum’s “‘Zwitter-ionen”’ 
theory of ampholytes might be possible. If this 
theory of ampholytes is correct, one would 
expect that gelatin, at its isoelectric point, should 
show a maximum dielectric constant; and the 
addition of either acid or base to isoelectric 
gelatin should cause a decrease. Again, the 
ultimate nature of a gel is a problem of con- 
siderable importance. What forces bind the water 
and maintain a solid condition even though the 
proportion of water is greater than 95 percent? 
No satisfactory answer has yet been made to this 
question. That the forces are electrostatic is often 
stated. A study of the dielectric constant of 
gelatin solutions under a variety of conditions 
offers a possible approach to this problem. If the 
water is bound, the dipoles are no longer free to 
orient and this part of the dielectric constant is 
decreased. 

Another possible application of dielectric 
constant measurements to proteins is the determi- 
nation of their apparent molecular weights. 

A review of the investigations on gelatin shows 
a contradiction in the results. Thus Drude! and 
Fiirth? found a decrease in the dielectric constant 
of gelatin solutions with increase in concen- 


1 Drude, Zeits. f. physik. Chemie 23, 267 (1897). 
? Firth, Ann. d. Physik 70, 63 (1923). 
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tration. These investigators, however, did not 
consider the effect of pH. Marie and Marinesco,' 
on the other hand, working apparently with 
more care in respect to pH of the gelatin and the 
manner of its preparation, found a rise of the 
dielectric constant with increase in concentration. 
The contradiction in these results has not been 
satisfactorily explained. It may be due to the 
difference in the frequency employed. Drude and 
Fiirth used a 76-centimeter wave while Marie and 
Marinesco used a six-meter wave. The possibility 
also exists that one or both sets of results do not 
represent true dielectric constants because of 
the uncertain influence of the conductivity of the 
solutions or other factors. 

While the present work was in progress, 
Marinesco* published a comprehensive paper on 
the proteins. Work has been done, also, by Fricke° 
and by May and Schaffer.® 


APPARATUS AND MATERIALS 


The methods used in this work are described 
in Papers I and II of this series.* A Drude 
method was used with the shorter waves and the 
resonance method with the longer waves. The 
range of wave-lengths available was from 3.5 to 
about 100 meters. 

’ Marie and Marinesco, J. chim. phys. 27, 455-70 (1930). 

4 Marinesco, J. chim. phys. 28, 51 (1931). 

5 Fricke, Kolloid Zeits. 56, 166 (1931). 

6 May and Schaffer, Zeits. f. Physik 73, 452 (1931). 

* J. Chem. Phys. 1, 836, 842 (1933). 
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Fic. 1. Effect of concentration on the dielectric constant 
of gelatin solutions as measured by the Drude method at 
various values of pH. 


The materials used were Baker’s commercial 
gelatin and Eastman’s isoelectric gelatin. Two 
samples of the latter were obtained, representing 
two different lots. The Baker gelatin contained 
about 1.5 percent ash, chiefly calcium salts. The 
Eastman product was practically ash-free. A 
glass electrode apparatus was used for measuring 
the pH. All solutions were prepared at 50°C. 
The pH was adjusted by adding NaOH or HCl in 
dilute solution to the solution of gelatin. 


DISCUSSION OF RESULTS 


A series of measurements was made at constant 
frequency (3.522m) with different concentrations 
of gelatin. The results obtained are shown by 
Fig. 1. The dielectric constants** are plotted as 
ordinates and the concentrations as abscissae. 
For comparison the results of Marinesco are 





** In place of the term, dielectric constant, the term, 
“square of the index of refraction,’’ should be substituted in 
this case, for the values given are uncorrected for any 
efiects of conductivity; and, in general, the n® term is 
larger than e. The term, dielectric constant, will be used for 
convenience, however, in the whole discussion. 
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Fic. 2. Effect of concentration on the dielectric constant 
of gelatin solutions at different values of pH. Results of 
Marinesco. 


given in Fig. 2. From these curves in Fig. 1 it 
may be observed that in general there is a 
decrease of the dielectric constant with increase 
in concentration of gelatin. Isoelectric gelatin, 
however, in the lower concentrations, gives a 
slight rise, as shown by curve 4. 

It may be observed, also, as is shown better in 
Figs. 3a and b, that the dielectric constant of 
gelatin appears to decrease on either side of the 
isoelectric point; similar results were obtained by 
Marinesco at low wave-lengths, as may be 
deduced from curves 3 and 4 of our Fig. 2 of his 
data. These results are in agreement with 
Bjerrum’s theory as outlined above. 

At the longer wave-lengths, it is found by 
comparison of the values obtained by the 
authors (Fig. 4) with those of Marinesco (Fig. 2, 
curves 1 and 2) that there is no agreement. 
Marinesco finds a sharp rise of the dielectric 
constant with concentration of gelatin to a 
maximum which appears at a concentration of 
1.5 percent in case of isoelectric gelatin, and at 
0.75 percent for a gelatin of pH 6.6. After the 
maximum is reached, the fall with concentration 
is rapid until the dielectric constants measured 
agree with those at the lower wave-lengths. 
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Fics. 3a and b. Effect of pH on the dielectric constant 
of gelatin solutions as measured by the Drude method at 
various concentrations. 
































95 ] | 
| 
9 | is 
10 
| eo 
K 
85 a ——_--s 
we: Az1 
€ 80; — | = 
Pf] 
15 | | | 
0 | 2 3 4 5 
Concentration, % 


Fic. 4. Effect of concentration and of wave-length on 
the dielectric constant of isoelectric gelatin solutions 
as determined by the resonance method. 


On the other hand, the values obtained by the 
present authors show no such sudden rises and 
attain no maximum within the concentration 
range studied, although the rate of increase of 
dielectric constant with concentration at the 
longer wave-lengths does fall off somewhat in the 
higher concentrations. Marinesco attributes the 
rise of dielectric constant to the effect of the 
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CASE 


dipoles of the gelatin molecules and the abrupt 
decrease to the union of dipoles to form nonpolar 
groups. Our results give no evidence of such a 
coupling; the slight flattening of the curves at the 
higher concentrations may properly be associated 
with the increased viscosity. 

As previously mentioned, dielectric constant 
studies should give some information concerning 
the association of gelatin molecules with mole- 
cules of the solvent. In an aqueous gelatin 
solution there are presumably four factors which 
may contribute to the total dielectric constant. 
They are (1) water molecules with the dipoles 
free, for which e=78.57 at 25°C; (2) water 
molecules the dipoles of which are fixed, «= 2.2; 
(3) gelatin molecules with dipoles not orienting 
in the field, e=5.4; (4) gelatin molecules, the 
dipoles of which are free. The gelatin aggregates 
are large and it should be possible to make the 
frequency so high that they cannot orient. If 
such were the case, the fourth factor above would 
disappear and the dielectric constant of the 
solution would be determined only by the other 
three. 

Following Marinesco, it may be shown that if 
N,, Nz and N; represent the number of grams of 
free water, of fixed water and of dissolved gelatin 
in 1 cc of solution; and d;, dz and d3 represent 
the corresponding densities, the following re- 
lationship exists: 








e—1/Ni No N53 ' . ; 
( +242) = Napit Nope Naps (1) 
e+2 


d, dz 4d; 
Let S=N,+Ne2=total g of water in 1 cc and 
consider d; = do, then 
e—1/S 


N 
—+—) =(S—Ne)pit Neopet Nsps. (2) 
e+2\d;_ dz 





pi, b2 and p; represent the specific polarizations 
of free water, of electrically saturated or fixed 
water and of solid gelatin, respectively. Since 
S/d,+N3/d3=1 (the unit of volume of the 
solution), then by substituting in (2) and 
transforming, 


N2=(SpitNsps—([(e—1)/(e +2) ])/(p1— b2)- (3) 


It should be possible by Eq. (3) to calculate the 
number of grams of water per cc, the dipoles of 
which are fixed, provided the p values are known 
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for fixed and free water and for gelatin with fixed 
dipoles. 

The necessary values were obtained and the 
amounts of fixed water were calculated from our 
own data at 3.522 meters, as well as from 
Marinesco’s data at 1.9 meters. These values are 
recorded in Table I. The first column gives the 


TABLE I. Fixation of water molecules by gelatin. 














Spec. e of 
vol. solu- ¢—! 
pH gelatin p; N3 S$ tion e+2 H 
6.6 .724 .4488 .005 .9964 76.4 .9617 131 
6.6 .724 .4488 .01 993 71.8 .9593 .270 
6.6 .724 .4488 .02 .9855 65.6 .9554 .110 


Baker's gelatin measurements by authors 


6.6 .790 .491 01015 .9920 77.1 .9621 —.2665 
6.6 .790 .491 .00507  .9961 77.9 .9621 —.1745 
6.6 .790 .491 .0203 .9842 76.3 .9617 —.289 
3.8 .790 491 .00507  .9961 74.5 .9608 .296 
3.8 .790 .491 .01015 .9920 73.0 .9600 0443 
Eastman’s gelatin measurements by authors 
6.6 .790 491 .01015 .9920 75.2 .9602 —.1328 
8.2. .790 .491 .01015 .9920 73.5 .9617 .0147 
3.8 .790 .491 .00446 .9965 76.3 .9617 .0663 
2.76 .790 .491 .00446 .9965 73.3 .9602 553 
Pi =0.9630 








pH of the solutions; the second, the specific 
volume of solid gelatin; the third, 3, or the 
specific polarization of solid gelatin; the fourth, 
the concentration of gelatin in grams per cc; the 
fifth, S, the total grams of water per cc; the 
sixth, e, the dielectric constant of the solution. 
The last column, H, contains the calculated 
values for the grams of water fixed per gram of 
gelatin, which may be simply calculated from N2 
of Eq. (3). In the use of-this formula »3; was 
calculated by using Fiirth’s value of 5.9 for the 
dielectric constant of solid gelatin. 

It is clear from our results, which give, in some 
cases, negative values for H, that this method 
cannot be used to determine the extent to which 
water is bound by gelatin. Although it is possible 
that the frequency used was not high enough 
completely to prevent the rotation of the gelatin 
molecules, a more likely interpretation of the 
results is that the water bound by the gelatin is 
still partially free to show its dielectric properties. 

The determination of molecular weights by 
means of dielectric constant measurements is 


based upon the equation given by Marinesco in 
the form 


M=3(d/C)(RT/nV), 


where 


M=molecular weight 

T =temperature; here 298°K 

n= viscosity of solvent; here 0.011 poises 

\= wave-length in region of dispersion 

R=gas constant; 83.2 x 10° ergs/degree 

V=specific volume of the material in solution; here 1.38 

cc/gm 
c=velocity of light, 3 10'° cm/sec. 


Since dispersion was found from 7 meters to 50 
meters, the molecular weights indicated range 
from 13,500 to 96,300. Gelatin appears, therefore, 
to possess molecules of different molecular 
weights. Although this method is yet unproved, 
the values obtained are reasonable at least. The 
use of the ordinary viscosity as the effective 
viscosity acting on the molecules is of course open 
to question. 

It is noteworthy that the dielectric constant of 
a gelatin solution is unaffected by the setting of 
the solution to give a solid. Solutions show a 
higher dielectric constant when first prepared and 
this gradually decreases, but the rate is not 
influenced by the setting, at least so far as may 
be observed. When a solution has visibly set, the 
dielectric constant continues to fall for some 
hours thereafter, although with a decreasing rate. 

These observations, while in agreement with 
those of Fricke,® are not in agreement with those 
of Marinesco.4 Marinesco assumes that the 
maximum of the dielectric constant-concen 
tration curve (Fig. 2, curve 1) determined at 6.5 
meters wave-length has a unique significance 
He attributes the rapid rise in dielectric constant 
of the solution to the dipole nature of the mole- 
cules of added gelatin. At the concentration 0.75 
percent, which is for the apex of the curve, he 
assumes the dipoles to combine and thus produce 
a rapid decrease in dielectric constant of the 
solution with increase in concentration. He 
further states in the case of solutions of pH 6.6 
that the apex of the curve represents the 
concentration at which gel formation is first 
possible. This explanation is not satisfactory 
because the same reasoning cannot be applied to 
isoelectric gelatin, since the maximum in this 
case is at about 1.5 percent concentration, while 
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it is well known that isoelectric gelatin will set 
before that concentration is reached and also that 
gelatin sets most readily at the isoelectric point 
and not least readily, as Marinesco’s interpre- 
tation would lead us to believe. 

The picture of the nature of gelatin solutions 
that may be drawn from our data is approxi- 
mately the conventional one. The gelatin mole- 
cules or particles in suspension are of graded size 
and of rather large molecular weights, ranging at 
least from about 5000 to 100,000. The gradation 
in size of molecule is apparently a continuous one, 


MALONE 
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although gradations in steps of the order 10,000 
are not excluded by the accuracy of the data 
obtained. This method does not exclude the 
possibility of other molecules or aggregates being 
present so long as they have no dipole moment. 
This situation is not likely in the case of proteins 
because of the large moments consistently 
obtained for the amino acids which are their 
decomposition products. It appears, also, that no 
coupling of the dipoles occurs and that gel 
formation is unaccompanied by any marked 
change in the dielectric properties. 
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Dielectric Constant Studies. IV. The Moments of Some Inorganic Compounds 


Mary G. MaAtone* AND A. L. FERGusoN, Chemical Laboratory, University of Michigan 
(Received December 11, 1933) 


The moments of some inorganic compounds have been 
determined, some in more than one solvent and the best 
values are believed to be: for SbCl;, 3.12; SbBrs, 2.47; SbIs3, 
1.58; AsI3, 0.96; AsBrz, 1.60; PI;, 0.0; AsF3, 2.65; and 
LiClO,, 7.84 Debye. The dielectric constant of solid AsF3; 
was measured and found to be 5.7. The moments of the 
molecules measured have been calculated from the 


electronegativity of the elements and compared with the 
measured values. A fair agreement is found with deviations 
greatest for molecules of highest polarizability and pos- 
sessing large moments. The electronegativity values were 
deduced from the measurements of the hydrogen com- 
pounds or from their periodic relationship. 





HE intimate relationship existing between 
the properties of a molecule and the 
electrical symmetry or electric moment has led 
to the belief that the latter gives considerable 
information concerning the nature and extent of 
the forces uniting atoms. Indeed so general has 
this belief become that the terms polar and 
nonpolar have often been confused with ionic and 
covalent. Atomic combinations between like 
atoms are generally known to be nonpolar and it 
is also generally accepted that true ionic bonds 
like those of the alkaline metal halides are 
strongly polar. A growing tendency may be noted 
to consider that there is no sharp division 
between bonds of the ionic and covalent type but, 
rather, that there is a general transition from one 
type to the other. 

It has been pointed out previously! that 
probably all covalent bonds are nonpolar and 
that the electric moment may be considered a 
measure of the deviation of the bonds from the 
covalent type. The measured values of the 
moments of simple molecules were so limited 
that the present study was undertaken in order 
to extend our knowledge in this direction and 
thus provide for a more adequate generalization 
in this field. 


* Rewritten from a thesis by Mary G. Malone in partial 
fulfillment of the requirements for the degree of Doctor of 
Philosophy at the University of Michigan. 

‘J. G. Malone, J. Chem. Phys. 1, 197 (1933). 


PREPARATION OF MATERIALS AND SOLUTIONS 


Benzene 


The thiophene-free material was stirred with 
sulphuric acid until it was no longer colored, 
washed with water and sodium hydroxide solution 
and stirred with water for some time. The 
benzene was then dried with calcium chloride and 
fractionated. The middle portion was fractionally 
frozen, redistilled from phosphorous pentoxide 
and dried over sodium. 


Carbon disulfide 


The c.p. quality was further purified by 
shaking with mercury until the mercury was no 
longer tarnished, fractionated, drie by refluxing 
with phosphorous pentoxide, again fractionated 
and redistilled just before using. No difference in 
capacity was noted in the different samples 
prepared if they were measured just after 
distilling. After standing for some time, however, 
a slight change was detected. 


Dioxane 


Eastman’s best quality was further purified by 
refluxing with sodium, fractionally distilling and 
standing over sodium. Different samples of 
dioxane showed different dielectric constants 
depending upon the length of time it was allowed 
to stand over freshly cut sodium. Constants for 
dioxane were redetermined every day it was used. 
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Carbon tetrachloride 


Eastman’s best quality was fractionally dis- 
tilled, dried over calcium chloride and redistilled. 


Antimony trichloride 

Baker’s c.p. grade was purified by fractional 
distillation under reduced pressure. The distillate 
was collected in the same vessel in which the 
solution was prepared so that no moisture was 
admitted. Benzene solutions so prepared showed 
no cloudiness. 


Antimony tribromide 

Coleman and Bell c.p. quality was purified in 
the same manner as the trichloride. The capacity 
of the carbon disulfide solutions changed with 
time, a noticeable change occurring in half an 
hour. To avoid this the solutions were made at 
25° and measurements taken immediately and at 
short intervals over a period of several minutes. 
The values obtained in this way did not change 
during the intervals. 


Antimony triiodide 

The Eimer and Amend c.p. quality was re- 
crystallized from carbon disulfide. The crystals 
were light in weight, bright red in color and 
every face shone with a luster. They did not 
tarnish when exposed to air for a short time. The 
carbon disulfide solution of antimony triiodide 
was clear and yellow in color and did not darken 
with time. The capacity measurements did not 
change with time. 


Arsenic tribromide 

The Eimer and Amend c.p. quality was 
purified by fractional distillation and collected in 
the manner described for the antimony chloride. 
The capacity measurements did not change with 
time. 


Arsenic triiodide 


The Coleman and Bell c.p. quality was 
purified by recrystallization from carbon di- 
sulfide. The crystals were red, shiny plates. They 
were dried from the carbon disulfide by dried 
carbon dioxide gas and handled in carbon 
dioxide during use. They were attacked by air, 
losing their luster. They formed an orange-to-red 
solution in carbon disulfide in contrast to the 
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yellow one produced by the antimony compound. 
The solution turned dark after several days. No 
change of capacity with time was noted during 
the time the solution was in the cell. 


Phosphorus triiodide 


It was prepared by refluxing, for some time, 
equivalent quantities of phosphorus and iodine 
dissolved in carbon disulfide. After the reaction 
was complete the solution was concentrated, 
cooled, and the red crystals were deposited. 
They were washed and recrystallized twice. The 
phosphorus triiodide fumed in the air and had to 
be handled in a dried carbon dioxide atmosphere. 
The carbon disulfide solutions were yellow in 
color and showed no tendency to darken when 
kept sealed from the air for three weeks. The 
capacity measurements were constant. 


Arsenic trifluoride 

It was prepared in a three-compartment dry 
glass apparatus which was evacuated. Equivalent 
weights of calcium fluoride and arsenious oxide 
were placed in the first compartment and a large 
excess of concentrated sulfuric acid added from a 
dropping funnel. The arsenic trifluoride formed 
was distilled into the second compartment, the 
third compartment being closed from it by a 
stopcock. The second compartment was closed 
from the first and the material kept until time 
for use. At that time a small first portion was 
distilled into the third compartment. It was 
removed by a vacuum pump and the middle 
portion distilled into this compartment for use. 
Solutions were made by adding the liquid from 
this compartment directly into the benzene in a 
dry nitrogen atmosphere. The glass apparatus 
and cell were only very slightly etched after use 
and the density bottle to a slightly greater 
extent. Some AsF; was kept in good condition in 
the glass apparatus for several months. 


Lithium perchlorate 


A very pure sample was obtained from 
Professor H. H. Willard. It was dried by fusing. 


Iodine monochloride 


Eimer and Amend c.p. quality was purified by 
a fractional recrystallization. 
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EXPERIMENTAL TABLE II. Moment of SbBrz in benzene at 25°. 


A heterodyne apparatus similar to that em- Po'=225 w= 3.28+0.1 


ployed by other workers? in this field was used for 
4 ‘ = C2 d , Pie P. Pra 
the solution measurements. For the solid arsenic - - 
trifluoride and iodine monochloride, a resonance 0.000 0.8374 = 2.276 
rn 1 0 ’ 0.00404 0.8858 2.351 27.60 50 


apparatus was used, described in a previous — 0.8899 2.369 27.82 
paper? from this laboratory. aairts ane — yi 


Two cells were used for the solution measure- 0.01270 0.9121 2.495 29.70 
ments, each consisting of five coaxial gold-plated 
brass cylinders. The plates of one of the cells 
were spaced with mica and those of the other 
with Bakelite. Leads were made to mercury cups. 
Different capacities could be obtained by using me d : Pw P» Pra Pr’ 
different combinations of the plates. The cells 
had maximum capacities of about 193 and 183 a eo] oe 21.69 180 50 130 
mmf, respectively. Small cells with platinum 9.003638 = =1.2745 2.692 21.85 170 120 


R 0.003959 1.2748 2.700 21.92 173 123 
plates were used with the resonance apparatus. 9997421 12888 2.757 2245 172 122 


The density bottle was similar to that de- 9.013659 1.3120 2.844 = 23.23 111 
scribed by Dr. Wooley.* The 33 cc bottle had a 2 
cc graduated pipette attached making it possible 
to take readings at different temperatures with 
one weighing. An error of 0.03 percent could occur 
= the volume reading. : : ss r d ’ P,. - Bas 

e apparatus was tested by measuring the — cS aan ROS 
moment of antimony trichloride in benzene. A 0.001603 1.2706 2.652 21.465 76.8 
‘ . : ...% 0.001945 1.2731 2.658 21.52 
value of 3.84 was obtained, which is in good 0.002140 12783 2600 21°815 
agreement with the values 3.6 by Werner,’ 4.11 0.002643 1.2797 2.665 21.56 
6 . 0.000 36° 1.2501 2.610 
ee and Engel® and 3.75 to 3.95 by 0003044 12778 2'662 21.70 
smith. 


In the Tables I to XI inclusive, c. represents 

















TABLE III. Moment of SbBr3 in carbon disulfide at 25°. 
Pop! = 225 w=247+0.1 

















TABLE IV. Moment of SbI; in carbon disulfide. 
Prp0/ = 52.5 w=1.58+0.1 

















TABLE V. Moment of AsBr; in carbon disulfide at 25°. 
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TABLE I. Moment of SbCl; in carbon disulfide at 25°. P20’ = 54.0 u=1.60+0.05 


© 








P20’ = 205 M =3.12+0.05 Ce d € Pr. P, Pra P,’ 








; ; 7 , "0.000 2.634 
- a Piz Pr Para 0.002648 2.651 21.478 0 37.6 
0.000 1.2600 0.005719 2.671 21.67 , 
0.00243 1.2659 21.85 238 44 0.008448 2.689 21.84 
0.00451 1.2685 22.20 ~3=—217 3 0.009160 2.697 21.91 
0.00707 1.2718 22.62 205 0.010454 2.707 21.90 
0.00839 1.2763 22.77 194 0.011953 2.713 22.07 
0.01281 1.2834 23.28 181 
0.01689 1.2930 24.14 187 
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TABLE VI. Moment of AsI; in carbon disulfide at 25°. 


*Smyth, Dielectric Constant and Molecular Structure, p. Pro’ = 19.0 u=0.96+0.1 
54, Chemical Catalogue Company, New York, 1931. 
‘J. G. Malone, A. L. Ferguson and L. O. Case, J. Chem. d € Pr P2 Priya Pr?’ 
Phys. 1, 836 (1933). 2.634 
* Wooley, J. Am. Chem. Soc. 55, 3518 (1933). 2.636 21.329 
* Werner, Zeits. f. Anorg. allgem. Chemie 181, 154 (1930). 2.648 21.42 


° Bergman and Engel, Zeits. f. physik. Chemie B13, 232 oo a 


(1931). 2.679 21.71 
‘Smith, Proc. Roy. Soc. (London) A136, 256 (1932). 











59 


“100 S777 
eS 100 we 
COMPU 











102 M. G. MALONE AND A. L. 


TABLE VII. Moment of PI; in carbon disulfide at 25°. 


u=0 








Ce € Pry P» Pra 


0.000 2.633 

0.001650 2.636 21.334 43.83 49.75 
0.002708 2.639 21.380 51.69 
0.004840 2.641 21.42 47.52 
0.005079 2.644 11.445 51.12 
0.007206 2.648 21.51 51.34 











TABLE VIII. Moment of AsF3 in benzene at 25°. 
Pr’ = 145.6 w=2.65+0.05 








C2 d P12 P, Pra —~P? 





0.000 0.8734 
0.01499 
0.02371 
0.02914 
0.03535 


0.8869 28.71 164.1 26.8 137.3 
0.8948 29.78 162.5 135.7 
0.8997 30.43 156.5 129.7 
0.9052 31.21 155.6 128.8 








TABLE IX. Moment of LiClO, in dioxane at 25°. 
Peco’ = 1285 u=7.84+0.05 








d « P12 P, Pe Ps 


1.0310 2.315 

1.0345 2.544 28.92 1048 24 
1.0362 2.620 29.92 990 
1.0370 2.672 30.40 994 
1.0310 2.325 

1.0318 2.381 26.88 1258 
1.0326 2.435 27.60 1172 
1.0336 2.498 28.37 1133 
1.0342 2.524 28.68 1106 











the mole fraction of the polar substance, ¢ the 
dielectric constant of the solution, d the density, 
P 2 the polarization of the mixture and P:2 the 
polarization of one solute. Pe is calculated from 
the following relationship, 


Pi3= [(e— 1)/(e+2) ]- (€:M,+¢2M?2)/d 


and 
P2= (P12—Pc1)/C2 


where ¢; and ¢c2, M; and M2, and P; and Pz» are 
mole fractions, molecular weights and polariza- 
tions of the components. 

The Pr,a values were obtained from the 
dielectric constants of the solid polar material 
obtained in case of SbCl;, SbBrs, AsBr;, AsI; and 
PI; from measurements by Schlundt.* The 


§ Schlundt, J. Phys. Chem. 5, 503 (1901); 8, 122 (1904); 
15, 381 (1911). 
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dielectric constant of solid AsF3 was measured by 
the use of the resonance apparatus and found to 
be 5.7. The Pg value of LiClO, was obtained from 
refractive index measurements. An Abbé re- 
fractometer was used; therefore this value 
represents the order of magnitude and not an 
exact measurement. For use in calculating the 
Px values the density of the solid material was 
obtained from the literature in all cases except 
PI; and AsF3;, which had to be measured. The 
density of PI; was found to be 3.89 at 25° and 
that of AsF;, 3.01 at the same temperature as the 
dielectric constant was measured. 

The P,.’ value represents P:—Pr+44. The 
values P.’ were obtained by extrapolating the 
P.2'—C2 curves to infinite dilution. The moment 
was obtained from the relation 


w=0.0127(P2.’'T)! Debye. 


An attempt is made to give an evaluation of 
the size of the known errors in giving the size of 
the moment. However, this is only an approxi- 
mation but given so that other authors may not 
make use of these values to a greater accuracy 
than the values justify. 

An attempt was made to measure the moment 
of HgCl. in benzene. Its insolubility made accu- 
rate measurements impossible, but the values 
indicated a zero moment, although a small 
moment would not be excluded. 

Considerable work was done in an attempt to 
measure ICI in CCl, solution. Just what percent 
of dissociation one has in such a case is not 
known. The red solution tarnished the gold plate 
of the cell during the set of measurements. The 
solution was kept in an inert atmosphere during 
use. The temperature and dielectric constant 
curve of solid ICI was studied by means of the 
resonance apparatus. A change of dielectric 
constant with temperature was marked, as is 
shown in Table X. The reason for such changes is 


TABLE X. The change of dielectric constant of solid ICI 
with temperature. 








Temp. Temp. Temp. 
—48° : — 28° , —15° 
—40° , — 25° : —12° 
—37° . — 20° \ —11° 
—34° ' — 18° , 
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not apparent. Somewhat similar changes in the 
dielectric constant of the hydrogen halides have 
been ascribed to the effect of rotational energy.* 


TABLE XI. Polarization of ICl in CCl, at 25°. 











Pro! = 60 

C2 € d Py. P, 
0.000 2.230 1.5837 
0.003423 2.239 1.5863 28.35 58.4 
0.005060 2.243 1.5880 28.38 55.9 
0.005126 2.244 1.5886 28.39 55.5 
0.006238 2.2465 1.5891 28.43 55.1 
0.008589 2.251 1.5912 28.47 54.7 








Strong electrolytic types of atoms are usually 
considered to possess high moments. The meas- 
urements given here for lithium perchlorate 
further substantiate this view. It may be that it 
is even a better case than the silver perchlorate 
measurement® so often quoted, for Gosh” was 
unable to show complex formation between 
perchlorates and dioxane. 

The theoretical treatment follows that of the 
work of Malone! previously referred to. The 
electronegativity scale has been extended on the 
basis of the periodic table, as is recorded in 
Table XII. The bonds not in parenthesis are 


TABLE XII. Electronegativity of non-metallic elements 
referred to H as zero. 











0.28 0.40 
N O F 
1.04 1.32 (1.72) 
0.69 
P S Cl 
0.36 0.63 1.03 
0.26 
As Se Br 
0.10 (0.37) 0.78 
0.40 
Sb Te I 
(—0.30) (—0.03) 0.38 








those taken from the experimental data on the 
hydrogen compounds. The bonds in parenthesis 
were obtained by the addition or subtraction of 
the differences in either the horizontal or vertical 
columns. For example, the value of F was 





* Profecsor Barker informs us that this is impossible in 
this case 
* Williams, Phys. Zeits. 29, 174 (1928). 
Gosh, J. Chem. Soc. 107, 159 (1915). 
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obtained by adding 0.69, the difference between 
N and P and between O and S, to the value of Cl 
to obtain 1.72 for F. 

In order to use such a table to calculate mo- 
ments it must be assumed that electronegativity 
is an additive property. Also, in the case of 
several vector moments present in a molecule, 
the angles present must be known. The reaction 
of one moment on another is a matter of some 
importance but its effect is difficult to evaluate. 
One assumes also that the angles remain constant 
as long as only two species of atoms are con- 
cerned. 

It is to be expected in highly polarizable atoms 
with high moments that considerable dis- 
crepancies would be encountered. But it may be 
assumed that the agreement will be the best in 
those cases where the conditions most nearly 
approximate the polarizability conditions preva- 
lent in the hydrides, upon which the table is 
based. This has proved to be the case and is to be 
observed, especially in the case of the arsenic 
fluoride. The agreement of the calculated and the 
observed value is very good. This value is one of 
great interest. All modifications of the scale given 
here indicate that fluorine is an element of greater 
electronegativity than any other. On this basis, 
the moment of HF should be greater than the 
moment of HCl. For the same reason the 
moment of AsF; should be greater than the 
moment of AsCl;, as was found to be the case. 
Some recent work by Kirkwood" has indicated a 
moment of about 0.8 for HF. Our results require 
a much higher value. A measurement of the 
moment of gaseous HF would be very valuable 
and is contemplated in this laboratory. 

Since complexes of benzene and antimony 
chloride and bromide are formed” it is of interest 
to have values for the moments of these materials 
in another solvent. Carbon disulfide was selected 
for the purpose, since there is no evidence 
available to indicate that complexes are formed 
with it. The measurements showed only slightly 
lower moments in each case. These lower values 
are listed in Table XIII. 

The antimony triiodide was measured only in 


carbon disulfide where its solubility is somewhat 


1! Kirkwood, Phys. Zeits. 33, 259 (1932). 
2 Smith and Davis, J. Am. Chem. Soc. 41, 411 (1922); 
Cremer and Duncan, J. Chem. Soc. 2750 (1930). 








































































































































































































































































TABLE XIII. Calculated and observed moments. 











Mol. Calc. Obs. Refr. Mol. Calc. Obs. 
SCl, 0.51 0.56 Smith SbCl; 2.23 3.12* 
PF; 2.01 SbBr;_ 1.81 2.47* 
PC; 0.98 0.8 Berg. SbI; 1.14 1.58* 

0.9 Smith 
PBr; 0.62 0.61 Berg. SbH; 0.50 
PI; 0.0 0.0* SbF; 3.41 
AsCl; 1.44 1.97 Berg. NCi; 0.0 
2.15 Smith 
AsBr;_ 1.05 1.66 | Berg. NBr; 0.36 
1.6 
AslI; 0.43 0.96* NI; 0.92 
AsF,; 2.53 2.65* NF; 0.97 
IBr 0.40 ICI 0.65 0:9-1.1* 








*Work in this laboratory. 

Smith, Proc. Roy. Soc. A138, 154 (1932); 136, 256 
(1932); Bergman and Engel, Zeits. f. physik. Chemie B13, 
232 (1931). 


greater than in benzene. The value found, 1.58, 
leads us to believe the value 0.40 ascribed to 
Williams by Werner" is too low. 

The measured moments of the three antimony 
compounds vary considerably from the calculated 
ones. It must be borne in mind, however, that 
these were calculated from the predictions in the 
periodic table and not on a measured hydride, as 
was true in the other cases. As was previously 
mentioned, molecules of high polarizability with 
high moments should be expected to deviate 
somewhat from the calculated values and it is 
borne out in these cases. In predicting a moment 
for SbF; from these data one would place it nearer 
the calculated value than the other antimony 
compounds. Its insolubility in nonpolar solvents 
so far available prohibited a measurement. 

Arsenic tribromide was measured in carbon 
disulfide in order to compare it with the measure- 
ments by Bergman and Engel in CCl,. It is 


138 Williams and Werner, Zeits. f. anorg. allgem. Chemie 
181, 154 (1930). 

44 Bergman and Engel, Zeits. f. physik. Chemie B13, 232 
(1931). 
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assumed that carbon tetrachloride and carbon 
disulfide are compounds of similar polarity. The 
agreement is all that could be desired. 

In the case. of PI; one would expect only a very 
small moment, as the moment of the bond P—I 
is negligible. This has been found to be the case, 
the measured moment could not be distinguished 
from zero. 

Any moment calculated from the data given 
for ICI (if one is justified in making such a 
calculation) is much higher than the one expected 
from the bond moments. ICI is another highly 
polarizable molecule with a reasonably high 
moment. 


CONCLUSION 





The evidence presented by the additional 
measured values of the dipole moments of 
inorganic compounds indicates the advantages 
and the limitations of attempting to treat the 
concept of moments on the basis of electronega- 
tivity. It is already indicated that the concepts 
which this method introduces are not wholly 
adequate to give quantitatively the moments in 
some cases. The method on the other hand has 
the merit of allowing generally recognized but 
poorly stated concepts to be expressed in terms 
that are capable of quantitative evaluation. A 
great mass of dipole moment data has been 
accumulated, chiefly on complex organic mole- 
cules, which if properly interpreted would, no 
doubt, be of great value. Such interpretations 
must be based on the fundamental relation of 
dipole moment to other properties of the mole- 
cule and these in the opinion of the authors can 
only be obtained by adequate study of the 
simpler cases. 

The authors gratefully acknowledge the as- 
sistance and interest of Dr. J. G. Malone. 
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This section will accept reports of new work, provided these 
are terse and contain few figures, and especially few halftone 
cuts. The Editorial Board will not hold itself responsible for 
opinions expressed by the correspondents. Contributions to 
this section must reach the office of the Managing Editor 
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not later than the 15th of the month preceding that of the 
issue in which the letter is to appear. No proof will be sent 
to the authors. The usual publication charge ($3.00 per 
page) will not be made and no reprints will be furnished free. 


The Diffusion of Hydrogen Through Metals: Fractionating the Hydrogen Isotopes 


In the summer of 1930 a striking example of the dif- 
fusion of hydrogen through metals was observed. The 
diffusion of the hydrogen took place under the following 
conditions. Two pieces of steel sheet metal, each about 
1/32’ (0.75 mm) thick, were welded together carefully at 
the four edges and immersed as cathode in an acid solution 
(mixture of sulfuric and hydrochloric acids). The tem- 
perature of the acid bath was about 70°C. The current 
density employed was about 22 amperes per dm*. At the 
end of two hours a very marked bulging of these sealed 
plates or iron “‘pillow’’ was observed. Thereupon the seal 
was broken and the gas collected over water and tested. 
The clear, sharp explosion that was obtained showed the 
presence of hydrogen. This experiment was repeated several 
times and always the same results were obtained. We 
decided to determine whether or not the hydrogen thus 
collected within the metal “‘pillow’’ was substantially the 
lighter of the two hydrogen isotopes. 

Preliminary qualitative experiments were carried out 
with tubes of iron, nickel and palladium. The walls of the 
iron tubes were about 0.8 mm thick, those of the nickel 
0.2 mm and 0.4 mm respectively, and those of the pal- 
ladium about 0.25 mm. The baths used were: a 25 percent 
(by volume) sulfuric acid solution, a 10 percent NaOH 
solution (by weight), and a 25 percent MgSO,-7H.O 
solution (by weight). The temperature varied from 85 to 
90°C. The current density employed varied from a little 
less than 11 amperes per dm? to over 120 amperes per dm’. 
The anode used was platinum. The gas passing through 
the walls of the tubes was collected over water. The pal- 
ladium tube gave the best performance. A 25 percent 
sulfuric acid solution was used and the gas gave a clear 
test for hydrogen. 

An apparatus was then set up whereby the hydrogen 
passing through the palladium tube wall could be “‘burned”’ 
and the resulting water collected. The electrolyte con- 
sisted of special water containing about } percent “heavy 
water,”’ to which was added fuming sulfuric acid, in the 
Proportion 335 cc special water, to 25 cc fuming sulfuric 
acid. During the electrolysis the volume of the solution was 


kept constant. The temperature was maintained at about 
85°C. A current density of 60 amperes per dm? was em- 
ployed. The voltage drop was 4.3 volts. A platinum anode 
was used. The length of the palladium tube immersed in 
the electrolyte was 25 cm. At the end of 28 hours, about 
2 cc of water had collected. Its specific gravity was found 
to be 1.00020. This would indicate that neither light nor 
heavy hydrogen alone passed through the walls of the 
palladium tube, but a mixture of both. 

An experiment was then arranged to determine the rate 
at which hydrogen passed through the walls of the pal- 
ladium tube used in the preceding experiment. The tube 
was carefully cleaned. A new solution was used consisting 
of distilled water 75 percent by volume; c.p. sulfuric acid 
25 percent by volume. The temperature of about 70°C was 
maintained for each run. The experiment was continuous, 
that is, with change of current density the apparatus and 
solution were left undisturbed. The volumes of gas were 
collected under as similar conditions as possible. The tem- 
perature of the water, over which the gas was collected, 
was 20.5°C. The barometric pressure was 766.4 mm. The 
data are given in Table I. 


TABLE I. 








Current Volume gas Volume gas 
density collected Time __ per minute 
amp./dm? Voltage cc min. cc 


100. 5 1000 48 20. 
50. ob 430 120 a 
20. 3.2 60 60 1. 











Another experiment was run to determine whether 
oxygen would diffuse through the walls of this same 
palladium tube. The same electrolyte was used as in the 
preceding experiment. A temperature of 40°C was main- 
tained. A current density of 100 amp. per dm? was used. 
The cathode was platinum. No gas was observed to come 
through the tube even after 2} hours. During the run, 
however, a dark violet-colored film was observed on the 
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surface of the tube. That film might have acted to prevent 
the diffusion of the oxygen through the walls of the tube. 

These preliminary experiments show that with the pal- 
ladium cathode fractionation of the two isotopes does 
occur, protium, the lighter isotope, diffusing through more 
readily. The experiments with the iron “pillow” indicated 
that we are dealing with extraordinary high pressures. On 
the basis of the overvoltages involved, a simple calculation 


Water containing approximately one part deuterium to 
200 parts of protium was electrolyzed with a current of 
one ampere and the hydrogen produced was passed over 
copper oxide at a temperature of 200°C. The rate of flow 
of the hydrogen was such that roughly 50 percent of it was 
converted into water and the remainder was forced over 
copper oxide at about 600°C, so that at the high temper- 
ature all the hydrogen reacted to form water. The two 
fractions of water, 10 ml each, were condensed separately, 
and distilled into pycnometers calibrated at 25°C. The 
sample of water formed at 200°C had a mass of 6.6 mg 
less than the 10 ml formed by the complete combustion. 


A recent letter by Coffin and Geddes? states that the 
first order coefficient for the decomposition of paraldehyde 
decreases by 50 percent when the pressure is increased 
from a few mm to 18 atm. They suggest as an explanation 
that the reactive degrees of freedom lie in the center of the 
molecule and are not excited directly by collisions; it is 
supposed that the probability of reaction per unit time 
will therefore increase when an activated molecule is 
undisturbed. In considering this mechanism it will be 
helpful to adopt the usual convention that the actual 
reaction may be mentally suppressed without otherwise 
changing the behavior of the molecule. States for which 
some small group of internal vibrations has at least a 
critical energy decompose spontaneously when the con- 
vention is relaxed; while it is maintained, they are present 
at the Maxwell-Boltzmann quota. Coffin and Geddes make 
the plausible-sounding assertion that the rate at which 
these states are produced is less at high pressures; to main- 
tain the same equilibrium, however, it would be necessary 
for the rate at which such states are destroyed to decrease 
to the same extent. It seems impossible to believe that the 
occurrence of collisions could prevent the central degrees 
of freedom from losing energy. A closer understanding of 
the fallacy is perhaps obtained by means of the following 
simple model; the molecule has three degrees of freedom, 
A, B and C; only A is directly affected by collisions, and 


We have previously pointed out! that the thermal dis- 
sociation of hydrogen iodide should be dependent on the 
concentration of deuterium in the hydrogen iodide. Ex- 
periments in this laboratory have verified these theoretical 
calculations. 
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Unimolecular Decomposition at High Pressures! 








Thermal Decomposition of Deuterium Iodide 





indicates that these pressures or the fugacities are of the 

order of a million atmospheres. 
CoLiIn G. FINK, 
Haro_p C. UREy, 
D. B. LAKE 

Columbia University, 
New York, New York, 
January 3, 1934. 





The experiment was repeated using two cells in series so 
that the rate of flow of hydrogen was twice as great and 
the surface of copper oxide was kept the same. (80 g of 
copper gauze, number 26 wire, was used in each reaction 
chamber.) This time the light fraction of 10 ml weighed 
5.0 mg less than the heavy fraction. The samples of water 
were interchanged in the pycnometers and identical mass 
differences were found. 
HERSCHEL HUNT 
Purdue University, 
West Lafayette, Indiana, 
January 11, 1934. 


A is connected with C only through the intermediary B; 
the probability that a quantum of energy will be trans- 
ferred from A to B, B to A, B to C or C to B is k sec... 
It is then readily shown that if A is given a quantum when 
B and C have none, the chance that C has the quantum 
at the expiration of time ¢ increases initially as (kt)*/2, and 
it would seem that such a model gives the Coffin-Geddes 
effect. This is not so, however; the collision will have no 
influence unless the quantum is at A, since only A par- 
ticipates in collisions; if the quantum is at A, the deac- 
tivation which occurs is perfectly balanced by activation 
of A in some other molecule. In general, it will be true that 
if collisions cannot directly excite some special vibration, 
they cannot directly quench it either. 
Louis S. Kasse.’ 
U. S. Bureau of Mines, 
Pittsburgh Experiment Station, 
Pittsburgh, Pennsylvania, 
January 12, 1934. 


1 Published by permission of the Director, U. S. Bureau 
of Mines. (Not subject to copyright.) 

2C. C. Coffin and A. L. Geddes, J. Chem. Phys. 2, 47 
(1934). 
3 Associate physical chemist. 


Hydrogen iodide was prepared from its elements and 
purified by three vacuum distillations. The reaction vessel 


1H. C. Urey and D. Rittenberg, J. Chem. Phys. 1, 137 
(1933). 
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was filled to about 60 cm pressure and sealed off. The reac- 
tion vessel was kept at the chosen temperature for a 
period about 15 times longer than would be necessary for 
ordinary hydrogen iodide to attain equilibrium. The reac- 
tion vessel was then pumped out through a break seal 
and the iodine and hydrogen iodide frozen out in traps with 
liquid air. The iodine and hydrogen iodide were dissolved 
in water and the solution analyzed for I, and H*. At no 
time were the gases in contact with anything but glass. 
The fraction decomposed, h, was calculated from the 
equation 


h=2(I2)/(2(1:J+(HI)), 


where [I.] is the moles of I; at equilibrium, and [HI] the 
moles of HI at equilibrium. The results are given in 
Table I. hg is the value Bodenstein? gives for ordinary 
hydrogen iodide. Ah is the difference in the fraction decom- 
posed of the sample rich in deuterium iodide and pure 
protium iodide 

The concentrations of the 13.5 and 40.0 percent samples 
were determined by W. Bleakney on his mass-spectro- 
graph. 

TABLE I. 








Ah calc. 
hp =0.20703 
0.00119 0.00175 
ha=0.22772 


0.00563 0.00602 
0.01153 0.01180 


h calc. Ah exp. 





0.20878 


0.20838 


T =468°C 
0.22724 
0.23287 
0.23877 


0.23326 
0.23904 
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The concentration of the 72.1 percent sample was 
determined by burning 500 cc of hydrogen iodide and 
determining the concentration of deuterium in the resultant 
water by the interferometric method which will be pub- 
lished shortly by Crist, Murphy and Urey. They assume 
that the molecular volume of deuterium oxide is the same 
as protium oxide. However, if the mol volume is that 
claimed by G. N. Lewis, the concentration of the sample 
will be about 76.5 percent. A change of 1 percent in the 
concentration of the deuterium will change Ah by 0.00013. 
This uncertainty will, however, not apply to the 13.5 and 
40.0 percent samples since the concentrations were deter- 
mined by the mass-spectrograph. It is to be noted that all 
our Ah’s are lower than those calculated. It is doubtful, 
however, if this is significant since we estimate that our 
experimental error is about 0.00040 in Ah. If we assume 
that our concentration of deuterium in the last sample is 
76.5 percent then the difference between the Ah exp. and 
Ah calc. will be about 0.00084. This seems to be definitely 
larger than our experimental error. 

While this experiment was conducted in the gaseous 
phase there is no reason to doubt that similar deviations 
from the classical value of the equilibrium constant will 
be found for reactions in condensed systems. 

D. RITTENBERG 
H. C. UREyY 
Department of Chemistry, 
Columbia University, 
New York, New York, 
January 15, 1934, 


2 M. Bodenstein, Zeits. f. physik. Chemie 29, 295 (1899). 
3 Bleakney, Phys. Rev. 34, 157 (1929); 39, 536 (1932). 


Heavy Water Inert, Due to Low Association 


Lewis and MacDonald! have reported the viscosities of 
heavy water over a considerable range of temperature. 
The very important question arises as to the effect of the 
extra proton in the nucleus of hydrogen upon the associ- 


TABLE I. 








50 100 
rhes rhes 


200 
rhes 





Hydrogen temperature constant 
(Bingham and Spooner) 

Oxygen temperature constant 

Observed temperature K for 
ordinary water 

Calculated temperature K for 
ordinary water 

Association of ordinary water 


45.9 
13.9 


269.7 


105.7 
2.55 


50.4 
13.6 


59.3 
13.0 
293.3 328.9 


114.4 
2.56 


131.6 
2.50 





Deuterium temperature con- 
stant 

Observed temperature K for 
deuterium oxide 

Calculated temperature K for 
deuterium oxide 

Association of deuterium oxide 


91.8 
277.8 


197.5 
1.41 


100.8 
301.7 


215.2 
1.40 


118.6 
(336.7) 


250.2 
(1.34) 








ation of water. Bingham and Spooner? have given the 
temperature constants for ordinary hydrogen and for 
oxygen at 50, 100 and 200 rhes, as reproduced in Table I. 
As fluidities have been found to be additive, the temper- 
ature required to give ordinary water a fluidity of 200 rhes 
would be 131.6°K, provided that the water remained unas- 
sociated. As shown in Table I, the observed value (328.9°K) 
is very much higher, leading to an association of 2.50. 
The fluidity of heavy water is quite different from the 
fluidity of ordinary water, so that it is necessary to conclude 
that the additional proton affects the fluidity of water and 
presumably therefore the temperature constant of fluidity. 
(See Fig. 1.) If this is so, the simplest assumption is that 
the effect is additive, so that if the constant of hydrogen is 
59.3 at 200 rhes the constant for deuterium is 118.6°, for 
the same fluidity. This would lead to a higher calculated 
temperature of 250.2° required to give deuterium oxide a 
fluidity of 200 rhes, again on the proviso that the water in 
question is non-associated. This is in accordance with the 
observed fact that heavy water is more viscous than 


1 Lewis and MacDonald, J. A. C. S. 55, 4730 (1933). 
2 Bingham and Spooner, J. Rheology 3, 221 (1932). 
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Fic. 1. Temperature-fluidity curve for deuterium oxide. 


ordinary water. But even though heavy water is more 
viscous, the same method of calculation indicates that 
deuterium oxide is very much less associated than is 
ordinary water, 1.34 as compared with 2.50. If the oxygen 
of ordinary water possesses ‘“‘residual valencies” as claimed 
by organic chemists, the additional protons in the deute- 
rium would be expected to lower the free energy of the oxide. 
Deuterium oxide is not only less associated than ordinary 
water but the temperature coefficient of association of 
deuterium oxide is less, as proved by the association at 50 
and 100 rhes. It has about the same association as such 





LETTERS TO THE EDITOR 


mildly associated substances as methyl formate or propi- 
onaldehyde. Deuterium oxide may therefore be regarded 
as rather inert chemically. The increased molecular weight 
would be expected to raise the freezing-point and boiling- 
point but the increase in the boiling-point given as 1.4° 
might be considerably greater were it not for the lower 
association. The decrease in surface tension, dielectric 
constant and refractive index are unanimous in indicating 
a lower association. The inert character of deuterium oxide 
is demonstrated by the fact that in it starch fails to react 
with enzyme, that seeds fail to germinate, that spirogyra 
fails to develop and that animal life perishes. It is assumed 
that deuterium oxide is not poisonous, since it is present 
in ordinary water, but it is quite clear that the life pro- 
cesses noted above may involve chemical reactions. 
Almost all aqueous solutions involve hydration. It is 
therefore not surprising that sodium chloride and barium 
chloride are relatively insoluble in heavy water, and it 
might be expected that the difference will prove to be less 
marked in the case of potassium chloride, silver nitrate 
and mercuric chloride, which dissolve with less evidence of 
hydration. 
EUGENE C. BINGHAM 
W. H. STEVENs, JR. 
Lafayette College, 
Easton, Pennsylvania, 
January 15, 1934. 
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